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Change  History 

•  Version  1.1:  First  publicly-available  version. 

•  Version  1.2:  Added  solutions  for  new  problems  (see  Problems,  Version  1.2  for  list). 
Corrected  solution  to  3.13-1. 
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Chapter  2 

Electric  and  Magnetic  Fields 


[m0002]  [1] 

2.2-1 


From  the  problem  statement,  V  =  1.5  V  and  d  —  30  // rn ,  so 


El  pa 


V 

~d 


50  kV/m 
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[mOOll]  [1] 

2.4-1 


From  the  problem  statement,  V  =  12  V,  er  =  6,  and  d  =  90  /mi,  so  the  electric  field 
intensity  is 

V 

|E|  «  —  ^  133  kV/rn 

LL 

Subsequently,  the  electric  flux  density  is 


D 


ereo  |E 


7.08  /iC/rn2 
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Chapter  3 
Transmission  Lines 
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[m0027]  [1] 

3.6-1 


(a)  The  expression  for  the  voltage  V(z)  traveling  in  the  +z  direction  contains  the  factor 
e~7Z.  The  propagation  constant  7  =  a  +  j/3,  where  a  and  (3  are  real- valued  constants. 
Therefore,  the  ratio  of  the  voltage  at  a  distance  /  from  some  other  point  on  the  transmission 
line  is: 


V{z  +  l) 
V(z) 


e—70H-0 

e~7Z 


=  e~71  =  e~ale~m 


The  magnitude  of  this  difference  is  just  the  hrst  factor;  i.e. ,  e  al.  We  also  know  that 


7  =  y/ {R'  +  juL')  [G'  +  jwC") 


At  100  MHz,  we  find  7  =  0.00850  +j3. 14468  m  T  Therefore,  a  =  0.00850  m  1,  and  the 
voltage  after  1  m  is 


(1  V)  exp  [—  (0.00850  m  x)  (1  m)] 


0.9915  V 


(b)  From  part  (a)  we  know  the  phase  of  this  difference  is  just  the  phase  of  the  factor  e  3/31. 
Since  /3  =  3.14468  rad/m,  the  phase  of  e~3®1  is  179.8°  for  l  =  1  m. 


(c)  For  a  radio  wave  in  free  space,  there  should  be  essentially  no  attenuation  over  1  m,  as 
long  as  this  1  m  span  is  located  far  from  the  transmitter.  This  is  because  free  space  propa¬ 
gation  contains  no  loss  mechanisms  analogous  to  R'  or  G'  in  the  transmission  line  model.  At 
/  =  100  MHz,  the  wavelength  of  the  radio  wave  is  A  =  c/f  =  3  m.  That  means  the  phase 
rotates  by  360°  in  3  m,  which  is  120°  in  1  m.  Note  that  the  wavelength  of  the  radio  wave 


is  significantly  longer  than  the  wavelength  of  the  signal  in  the  transmission  line. 


[m0027]  [2] 

3.6-2 


The  question  is  whether 

V(z)  =  V0+e~7Z  +  V0~e+7Z 
is  a  solution  to  the  TEM  transmission  line  wave  equation 

<92  ~  ~ 

—  V(z)-7aV(z)  =  0 

where  Vq+,  and  7  are  complex- valued  constants.  To  determine  this,  we  substitute  the 
candidate  solution  into  the  equation  and  determine  if  equality  holds.  Taking  the  first  deriva¬ 
tive  of  the  candidate  solution: 

^V(z)  =  -7U+e-^  +  yVo-e+r* 

Repeating  to  get  the  second  derivative: 

d2  ~ 

g^V(z)  =  +7V.V”  +  72r0“e+7a 
Now  making  the  substitutions  into  the  left  side  of  the  wave  equation: 


l+'fV+e-'*  +  72^0“e+7"]  -  72  |%+e~7Z  +  V0~e+^z] 

=  -^V+e-7*  +  72R0-e+72  -  72Vr0+e-72  -  72R0'e+72 

=  0 

which  is  the  the  right  hand  side  of  the  wave  equation,  as  expected. 
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[m0052]  [1] 

3.7-1 


It  is  true  that  the  real  part  of  the  characteristic  impedance  must  be  positive. 


Here’s  a  mathematical  argument:  Recall: 


Zn  = 


R!  +  jojL' 
G'  +  juC' 


Also  note  that  R! ,  L',  G'  and  C  must  all  be  positive  or  zero.  Therefore,  the  numerator 
and  denominator  of  the  above  expression,  before  taking  the  square  root,  must  both  have 
phase  in  the  range  0  to  +7t/2.  This  means  the  numerator  divided  by  the  denominator,  again 
before  taking  the  square  root,  must  have  have  phase  in  the  range  — 7t/2  to  +7t/2.  Taking 
the  principal  square  root  reduces  the  phase  by  a  factor  of  two,  the  phase  of  Z0  must  be  in 
the  range  — 7t/4  to  +7t/4.  Subsequently,  the  real  part  of  Z0  must  be  positive. 


Here’s  a  physical  argument:  A  positive- valued  real- valued  component  of  an  impedance  rep¬ 
resents  the  dissipation  of  power  (e.g.,  resistance)  or  transfer  of  power  out  of  a  system  (e.g., 
to  a  load).  Conversely,  a  negative- valued  real- valued  component  of  an  impedance  represents 
the  creation  of  power  or  the  introduction  of  power  into  a  system;  in  other  words,  an  active 
device.  Since  the  concept  of  characteristic  impedance  applies  to  transmission  lines,  and  since 
transmission  lines  are  passive  devices,  the  real  component  of  the  characteristic  impedance 
must  be  positive. 
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[m0080]  [1] 

3.8-1 


The  physical  current: 

i(z,  t )  =  (2  A)  sin((3  rad/s)t  +  (4  rad /m)z  +  5  rad) 

=  (2  A)  cos((3  rad/s)t  +  (4  rad/m)z  +  5  rad  —  n/2) 


so 


I(Z)  =  [(2  A)  ej^{b~n/2)  rad)l  ei(4  rad/m>* 
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[m0080]  [2] 

3.8-2 


Converting  to  time-domain  representation: 

v(x,t)  =  Re  {v(x)ejujt}  =  Re  {V0e+j^xejujt} 

The  problems  statement  implies  Vo  is  complex-valued.  To  accomodate  this,  we  define  the 
magnitude  and  phase  of  Vo  as  follows: 

V0  =  |  Vb  |  e^'3 


Then: 

v(x,t)  =  Re  { |  Vo  |  ej7T/3e+jPxeju}tj  =  \V0\  Re  {e^+^+7r/3)} 
Finally,  using  the  identity  e^6  =  cos  9  +  j  sin  0,  we  obtain 


v(x,  t ) 


|  Vq  |  cos  (cot  +  j3x  +  7t/3) 


This  wave  is  traveling  in  the  —x  direction. 
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[m0080]  [3] 

3.8-3 


The  form  given  in  the  problem  statement  is  a  phasor,  describing  a  wave  traving  in  the 
+0  direction.  To  obtain  the  time  domain  form: 

v(4>,  t)  =  Re  {V0e~j^ejuit}  —  |Vo|  cos  (cut  —  /3(f>  +  /) 

where  /,  the  phase  of  Vo,  is  not  given.  Note 

27 r 

/3  =  —  =  62.832  rad/m 
A 

Furthermore,  we  are  told  that  u(0  =  A/4,  t  =  0)  is  a  maximum,  so: 

I  Vo  |  cos  =  |  Vb  |  cos  (-|  + 

is  a  maximum,  which  means  =  +7t/2.  Therefore: 


v(<j>,t) 


|  Vo  |  cos  (ut  —  [62.832  rad/m]  </>  +  —  j 


The  problem  statement  does  not  provide  sufficient  information  to  determine  |Vq|  or  cu. 
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[m0080]  [4] 

3.8-4 


Since  one  end  of  the  transmission  line  lies  at  infinity,  we  expect  a  wave  traveling  in  the  +z 
direction  only.  (Note  for  future  reference:  The  same  effect  can  be  achieved  for  a  finite-length 


line  by  perfectly  matching  the  transmission  line  at  the  end  opposite  the  voltage  source). 


The  general  form  for  the  physical  (real-valued)  voltage  wave  is 


v+(z,  t)  =  Vq  e  az  cos  (cot  —  /3z  +  ip) 


(3.1) 


Examining  the  problem  statement,  we  determine: 

|  Vq  |  =  2  mV 

v+(z  =  0,  t  =  0)  =  —  2  mV 

|V+|  e-“'(10  =  1  mV 

/  =  15  MHz  (frequency  of  the  source) 

Wavelength  in  the  line  A  =  0.4Ao  where  Ao  is  wavelength  in  free  space. 

Quantities  remaining  to  be  determined  in  Equation  13.  II  are:  attenuation  constant  a,  angular 
frequency  co,  phase  propagation  constant  /3,  and  wave  phase  reference  ip. 

The  attenuation  constant  may  be  determined  as  follows: 


Vq+|  e““'(1°  m)  e-“'(1°  m)  1  mV  1 


(3.2) 


V+|e-“-(°m)  i  2  mV  2 


Therefore 


Solving  for  ck,  we  obtain  a  =  0.0693  m  1 


The  angular  frequency  is  simply  2irf  =  94.2  Mrad/s  =  co 


The  phase  propagation  constant  is  determined  as  follows: 


(3.3) 


The  wave  phase  reference  is  determined  as  follows: 

v +(z  =  0,  t  =  0)  =  |  Vq  |  •  1  •  cos  (0  —  0  +  i/j)  =  —  2  mV 


(3.4) 


Solving  for  ip,  we  find  ip  —  n 


The  boxed  quantities  above  comprise  the  complete  solution  to  part  (a). 
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The  solution  to  part  (b)  -  the  phasor  representation  -  is  simply: 


In  this  case,  we  obtain: 


since  eJ7r  =  —  1. 


(3.5) 

(3.6) 
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[m0083]  [1] 

3.9-1 


From  the  problem  statement,  Z0  =  72  fi,  L'  =  0.5  /xH/m,  /  =  80  MHz,  and  the  low-loss 
approximations  apply.  Using  the  low-loss  approximation  Zq  L' /C'\ 


C 


L' 

Z$ 


=  96.4  pF /m 


Subsequently,  the  phase  velocity  is 


1 


vp  ~ 


VI7cJ 

and  the  phase  propagation  constant  is 


1.44  x  108  m/s 


j3  ~  ojV L'C  =  27 \f\fl7C'  =  3.49  rad/m 
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[m0143]  [1] 

3.10-1 


The  characteristic  impedance  Zq  of  coaxial  cable,  assuming  the  low-loss  assumptions  ap¬ 
ply,  is 


where  er  is  the  relative  permittivity  of  the  spacer  material,  and  a  and  b  are  the  radii  of  the 
inner  and  outer  conductors,  respectively.  Air  has  er  ~  1  and  is  lossless  to  a  very  good  ap¬ 
proximation.  We  are  also  told  the  resistance  of  the  inner  and  outer  conductors  is  negligible. 
Therefore,  the  low-loss  assumptions  apply,  and  we  are  justified  in  using  the  above  expression. 

One  way  to  reduce  Z0  from  90  O  to  62  O  is  to  replace  the  air  spacer  with  a  material  spacer 
having 


Thus,  one  solution  is  to  replace  air  with  a  low-loss  material  having  er  =  2.11.  Another  way 
is  to  reduce  b/a.  For  the  90  O  cable,  we  determine  that 


To  reduce  Z0  to  62  0,  we  require 


Thus,  a  second  solution  is  to  keep  air  as  the  spacer  material  but  reduce  b/a  to  2.81. 
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[m0143]  [2] 

3.10-2 


Under  low-loss  conditions,  the  characteristic  impedance  Z$  of  a  coaxial  cable  is  given  by 


(3.7) 


where  er  is  the  relative  permittivity  of  the  spacer  material,  6  is  the  radius  of  the  outer 
conductor,  and  a  is  the  radius  of  the  inner  conductor.  Since  geometry  may  not  be  changed, 
In  (6/a)  may  not  change.  The  only  free  parameter  remaining  is  er.  Z0  is  maximized  by 
minimizing  er.  Since  the  minimum  practical  value  of  er  is  «  1  (i.e. ,  air,  or  a  vaccuum),  the 
optimal  new  value  of  er  is  1.  This  increases  Zq  as  follows: 


(75  ft) 


\/T25 


112.5  ft 


(3.8) 
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[m0084]  [1] 

3.12-1 


The  voltage  reflection  coefficient  is 


Zl  —  Z0 
Zl  +  Z0 


500  ft  -  75  ft 
500  ft  +  75  ft 


=  +0.739 


Therefore,  the  peak  voltage  of  the  reflected  wave  at  the  antenna  input  is 


(0.739)  (30  V)  +  22.2  V 


The  line  is  lossless,  so  there  is  no  attenuation  of  the  reflected  wave  along  the  return  trip  from 
antenna  to  transmitter.  Therefore,  the  peak  voltage  of  the  reflected  wave  at  the  output  of 
the  transmitter  is  22.2  V. 
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[m0084]  [2] 

3.12-2 


From  the  problem  statement,  V^+  has  magnitude  7  mV  and  phase  180°,  so  V^~  =  —7  mV 
Also  from  the  problem  statement  Zq  =  60  and  Zl  =  20  Therefore, 


Zl  ~  Z0 
Zl  +  Z0 


-0.5 


Subsequently,  V0  =  rV0+  =  +3.5  mV.  Thus,  the  magnitude  of  the  reflected  wave  is  3.5  mV 
and  the  phase  is  07 
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[m0084]  [3] 

3.12-3 


33  Q  -  140  n 


The  voltage  reflection  coefficient  is 

p  _  ZL  -  Z0  _  _ _ 

“  ZL  +  Z0  -  33  Q  +  140  n  ~ 

Therefore,  the  magnitude  of  the  reflected  voltage  wave  is 


|T  (3  V)|  =  1 1.86  V 

and  the  phase  of  the  reflected  voltage  wave  is 


170°  +  180° 


-10c 


0.6185 
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[m0084]  [4] 

3.12-4 


In  general,  the  voltage  reflection  coefficient  T  for  a  load  impedance  ZL  connected  to  a  trans¬ 
mission  line  having  characteristic  impedance  Z(l  is 


Zl  ~  Z0 
Zl  +  Z0 


Solving  for  ZL,  we  obtain 


ZL  = 


i  +  r 
l-r 


For  r  =  0,  the  formula  gives  Z0l  as  expected. 
For  r  =  +1,  the  formula  — >  oo,  as  expected. 
For  r  =  — 1,  the  formula  gives  0,  as  expected. 
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[m0086]  [1] 

3.13-1 


(a)  The  current  at  a  voltage  maximum  is  zero,  (b)  The  voltage  at  the  short  circuit  ter¬ 
mination  is  zero.  The  distance  between  voltage  extrema  is  A/4,  so  A/4  =  8  cm.  The  distance 
between  voltage  maxima  is  A/2  =  16  cm.  Therefore,  the  distance  between  the  short  circuit 
and  the  second  voltage  maximum  is  8  +  16  =  24  cm. 
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[m0081]  [1] 

3.14-1 


First  note 

So  in  this  case 

Also  note: 


where  in  this  case  Zq 
we  find: 


SWR  —  1 
SWR+  1 


in  < 


1.2 


1.2  +  1 


+  0.091 


p  _  ZL  -  Zq 

Zl  +  Zq 

50  hi  and  Zjj  is  the  input  impedance  of  the  amplifier.  Solving  for  Zl 


Zl  = 


Z o 


i  +  r 
l-r 


Since  the  imaginary  component  of  Z0  is  zero,  and  since  the  imaginary  component  of  Zl  is 
negligible,  T  must  be  real-valued.  Therefore,  —0.091  <  Y  <  +0.091  and 


4i.7  n<zL<  6o.o  n 
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[m0081]  [2] 

3.14-2 


From  the  problem  statement, 
tion  coefficient  is 


Zq  =  72  and  ZL  =  60 


r 


Zl  ~  Zq 
Zl  +  Zq 


-0.091 


Therefore,  the  voltage  reflec- 


and  the  standing  wave  ratio  is 


SWR 


r 

1  2 

|i'| 
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[m0081]  [3] 

3.14-3 


From  the  problem  statement,  Z0  =  50  and  Zl  =  20  —  j35  Therefore,  the  voltage 
reflection  coefficient  is 

Zl  —  Z0 


r  = 


Zl  +  Zq 


-0.143-  j'0.571 


and  the  standing  wave  ratio  is 


SWR  = 


i  +  |r| 
i  -  in 


=  3.87 
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[m0087]  [1] 

3.15-1 


The  input  impedance  of  a  lossless  line  is  periodic  in  length,  with  period  A/2.  Therefore, 
the  line  is  exactly  3  periods  long,  which  means  the  input  impedance  is  equal  to  the  load 
impedance  72  +  j'42 
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[m0087]  [2] 

3.15-2 


From  the  problem  statement:  Zq  =  50  fl,  Zl  =  Rdut  —  10  O,  and  l  =  10  cm.  Also, 
the  wavelength  in  the  transmission  line  A  =  0.6Ao,  where  Ao  is  the  free-space  wavelength.  As 
always,  (3  =  2tt /A  and  \o  =  c/ f  where  c  is  the  speed  of  light  in  free  space. 


Here’s  the  result  (see  end  of  this  solution  for  source  code): 


■  Re(Z) 
•  lm(Z) 


-200  L 
106 


108 

Freq  [Hz] 


The  answers  to  parts  (b)  and  (c)  depend  on  one’s  interpretation  of  “significance.”  Two 
reasonable  interpretations  are  (1)  a  qualitative  judgment  based  on  when  the  curves  seem  to 
clearly  diverge  from  the  nominal  (DC,  or  equivalently  l  =  0)  value  and  (2)  a  quantitative 
judgment  based  on  when  the  real  part  is  in  error  by  more  than  5%  (or  some  other  percentage) 
and  the  imaginary  part  is  in  error  by  more  than  5%  of  the  real  part.  Here  are  the  results 
using  both  criteria: 


Nominal  (l  =  0) 

“Qualitative”  >  5%  error 

RealjZ} 

io  n 

~  100  MHz  =  6.4  MHz 

Imag{Z} 

o  n 

~  10  MHz  =  3.0  MHz 

In  both  cases  it  is  clear  that  error  in  the  imaginary  part  is  significantly  degraded  at 
a  lower  frequency  than  the  error  in  the  real  part,  and  that  both  are  exhibiting 
large  errors  at  frequencies  greater  than  ~  10  MHz. 

Here’s  source  code  in  Octave  (should  also  work  in  MATLAB): 
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clear  all; 
close  all; 

ZL  =  10.0;  7,  [ohm]  R_DUT 

Z0  =  50.0;  l  [ohm]  characteristic  impedance 
1  =  0.1;  7.  [m]  length  of  line 

n=0;  7.  counting  points 

for  logf=6:  .001 : 10,  7.  incrementing  frequency  in  log  scale  from  10“6  to  10“10  Hz 
n=n+l ; 

f (n)  =  10.“logf;  7.  [Hz]  frequency 

lambdaO  =  (3.0e+8)/f (n) ;  %  [m]  free  space  wavelength 

lambda  =  0.6*lambda0;  l  [m]  wavelength  in  line 

b  =  2*pi/lambda;  7,  [rad/m]  beta  =  phase  propagation  constant  in  cable 

Z(n)  =  Z0*(ZL+j*Z0*tan(b*l))/(Z0+j*ZL*tan(b*l)) ; 

end 

semilogx(f ,real(Z) , ’b-’) ;  hold  on; 
semilogx(f ,imag(Z) , ’r-’) ;  hold  off; 
legend( ’Re(Z) ’ , ’lm(Z) ’) ; 
grid  on; 

xlabel ( ’ Freq  [Hz] ’ ) ; 
ylabel ( ’ Z  [ohm] ’ ) ; 

[f’  real(Z)’  imag(Z)  ’  (real(Z) ’-10)/10  imag(Z)’/10]  7.  used  to  answer  parts  (b)  and  (c) 
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[m0087]  [3] 

3.15-3 


From  the  problem  statement:  Z0  =  50  and  Zl  =  25  +  j 25 

(a)  Voltage  reflection  coefficient: 

p  _  ZL  +  Z0 
Zl  ~  Z0 

(b)  The  input  impedance  may  be  calculated  using 

1  +  Te~j2/31 

Zin  =  Z°l  -  pe— 42/3Z  (3.10) 

where  (31  =  (2n/\)l  =  2n{l/\).  The  requested  plot  is  shown  below.  In  this  figure,  “o” 
indicates  l  =  0  and  “x”  indicates  l  =  0.45A.  (See  end  of  this  solution  for  source  code.) 


-0.2+  jOA 


(3.9) 
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(c)  Here  are  the  lengths  for  which  the  input  impedance  is  completely  real-valued: 


l  =  0.162A 

=  130.9  n 

l  =  0.412A 

— >■ 

=  i9.i  n 

Here’s  source  code  in  Octave  (should  also  work  in  MATLAB): 

cleax  all; 
close  all; 

ZL  =  25.0+j*25.0;  7.  [ohm] 

ZO  =  50.0;  '/,  [ohm]  characteristic  impedance 

Gamma  =  (ZL-ZO) ./(ZL+ZO)  l  voltage  reflection  coefficient 
n=0;  */,  counting  points 

for  1=0:  .001:0.45,  l  [lambda]  incrementing  length  from  0  to  almost  lambda/2 
n=n+l ; 

bl  =  2*pi*l;  7.  [rad]  electrical  length 

Z(n)  =  Z0*(l+Gamma*exp(-j*2*bl))/(l-Gamma*exp(-j*2*bl)) ; 

end 

hi  =  plot (real (Z) ,imag(Z)) ; 
axis ("equal") ; 
grid  on; 

xlabel( ’Re [Z]  [ohm]’); 
ylabel ( ’ Im  [Z]  [ohm] ’ ) ; 

1=0.00;  7,  [lambda] 

bl  =  2*pi*l;  l  [rad]  electrical  length 

Zp  =  Z0*(l+Gamma*exp(-j*2*bl))/(l-Gamma*exp(-j*2*bl)) ; 

hold  on;  h2  =  plot (real (Zp) ,imag(Zp) , ’ro’) ;  hold  off; 

1=0.45;  7.  [lambda] 

bl  =  2*pi*l;  7.  [rad]  electrical  length 

Zp  =  Z0*(l+Gamma*exp(-j*2*bl))/(l-Gamma*exp(-j*2*bl)) ; 

hold  on;  h3  =  plot (real (Zp) ,imag(Zp) , ’rx’) ;  hold  off; 
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[m0088]  [1] 

3.16-1 


In  this  case,  the  input  impedance  is 


Z stub  jZo  COt  [31 


where  Z0  =  75  fl,  l  =  13  cm,  and 


13  —  —  —  2nf 


vp  0.55c 


where  /  =  900  MHz.  Therefore,  (3  =  34.3  rad/m,  [31  =  4.45  rad,  and  Zstub  =  —j  19.7  H 
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[m0088]  [2] 

3.16-2 


From  the  problem  statement:  Z0  =  75  f  =  1.5  GHz,  Zin  =  +j'300  fl  is  desired,  and 
vp  =  0.6c.  Note  that  for  a  short  circuit,  in  this  case: 

Zin  =  +jZ{)  tan  fjl  =  +j300  fl 


so 

Note 

so  l  =  2.53  cm 


f3l  ^  1.3258  rad 

(3  =  —  =  =  52.36  rad/m 

vp  0.6 c 
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[m0088]  [3] 

3.16-3 


From  the  problem  statement:  /  =  5.8  GHz,  Z$  =  50  vp  =  0.7 c,  and  the  capacitor  to 
be  replaced  has  value  C  =  83  pF.  Therefore  the  desired  impedance  is 


^c  = 


3 


27 TfC 


^  -j'0.3306  n 


We  choose  an 


open-circuited  line,  as  this  yields  the  negative  reactance  for  the  shortest 


possible  lengths.  The  input  impedance  of  an  open-circuited  line  is 

Zin  =  —jZ0  cot  (31 

Setting  this  equal  to  Zc  and  solving  for  (31: 

PI  =  cot"1  Q  =  1-5642  rad 


Note 


(3=-  = 


-j  50  n 

c o  2i xf 


vp  0.7c 


=  173.54  rad/m 


so 


l  =  9.01  mm  . 
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[m0145]  [1] 

3.17-1 


(a)  For  a  bandpass  response  centered  at  200  MHz,  you  want  the  input  impedance  into 
the  stub,  which  is  attached  in  parallel  to  the  line,  to  be  an  open  circuit  at  200  MHz.  This  is 
accomplished  using  a  short-circuited  stub  which  is  one  quarter  wavelength  long  at  200  MHz. 
A  wavelength  in  the  transmission  line  is 


A  = 


0.67  (3  x  108  m/s) 
200  MHz 


=  1.005  m 


(3.11) 


so  the  stub  length  is  25.12  cm  . 


(b)  See  below: 
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[m0145]  [2] 

3.17-2 


The  smallest  length  for  which  the  imaginary  component  of  the  impedance  of  an  open- 
circnited  stub  is  positive  is  slightly  greater  than  A/4.  The  imaginary  part  of  the  impedance 
remains  positive  until  the  length  is  slightly  less  than  A/2.  In  this  transmission  line, 


0  = 


U ! 
Vp 


27 r/ 
0.7c 


89.8  rad/m 


where  /  = 


and  so  the 


]  GHz.  Therefore, 

\  2n  7 

A  =  —  =  7  cm 

r 

smallest  contiguous  range  of  transmission 


1.75  cm  <  l  <  3.5  cm 
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[m0145]  [3] 

3.17-3 


For  zero  response  centered  at  fc=  1.3  GHz  you  want  the  input  impedance  into  the  stub, 
which  is  attached  in  parallel  to  the  line,  to  be  a  short  circuit  at  /  =  fc.  This  is  accomplished 
using  an  open-circuited  stub  that  is  one  quarter  wavelength  long  at  /  =  fc.  A  wavelength 


in  the  transmission  line  is 


,  0.6  (3  x  108  m/s) 

A  =  1.3  GHz  =  13  84  Cm 


(3.12) 


so  the  stub  length  is  3.46  cm  .  The  characteristic  impedance  is  irrelevant. 
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[m0091]  [1] 

3.19-1 


From  the  problem  statement,  we  see  that  the  design  will  consist  of  a  quarter- wave  matching 
section  followed  by  a  line  having  a  characteristic  impedance  Zq2  =  300  hi  -  i.e. ,  equal  to  the 
load  impedance  -  and  the  total  length  will  be  l  =  5  cm.  The  characteristic  impedance  of 
the  quarter-wave  section  must  be  Zm  =  \fZsZ02,  where  Zs  is  the  source  output  impedance; 
thus,  we  have  Z01  =  122.5  12.  The  length  l\  of  the  quarter-wave  section  is  A/4,  where  A  is 
the  wavelength  in  the  transmission  line.  For  FR4,  we  have 

A  =  =  c//  =  (3  x  10°  m/s)  /  (1.5  GHz)  =  ^  ^ 

yi(Cr  +  l)  /i  (4.5  +  1) 

so  l\  =  3.01  cm  and  l2  —  l  —  h  —  1.99  cm. 


What’s  left  to  figure  out  is  the  width  w  of  the  microstrip  lines,  which  determines  the  charac¬ 
teristic  impedance  since  h  =  1.6  mm  and  er  =  4.5  are  already  set.  We  know  that  h/w  =  1/2 
gives  a  characteristic  impedance  of  5011  for  FR4,  so  the  width  of  a  5011  line  is  2 h  =  3.2  mm. 
To  get  the  higher  characteristic  impedance  ^01  =  122.511,  W\  will  have  be  smaller  than 
3.2  mm.  An  approximate  but  reasonable  solution  is  simply  to  assume  the  characteristic 
impedance  scales  with  w  in  the  same  way  (i.e.,  linearly)  as  it  does  in  the  “wide”  ( h/w  <C  1) 
case,  so 

w/i  ~  (3.2  mm) - =  1.3  mm 

v  ;  122.5  12 


and 


w2  & 


(3.2  mm) 


50  12 
300  11 


0.5  mm 


You  could  also  use  the  Wheeler  (1977)  formula  or  some  other  equation  or  reference;  how¬ 
ever,  the  increased  accuracy  is  typically  irrelevant  in  practice  due  to  issues  such  as  the  large 
variation  in  er  due  to  manufacturing  issues.  So,  while  it’s  not  wrong  to  take  that  approach, 
it’s  usually  not  worth  the  effort  if  you  are  able  to  instead  “scale”  from  a  known  design  as  we 
have  done  above. 


So,  your  sketch  should  show  the  source,  followed  by  3.01  cm  of  line  which  is  1.3  mm  wide, 
followed  by  1.09  cm  of  line  which  is  0.5  mm  wide,  followed  by  the  load,  as  shown  in  the  figure. 


n 

A 

[So^ 


r 


Off 
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[m0091]  [2] 

3.19-2 


From  the  problem  statement:  Z l  =  200  l  =  A/4,  and  Z0  =  100  £2.  Since  this 
quarter- wave  line, 

L 


7.  _  _ 0  _ 

^in  ry 


50  n 


is  a 
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[m0091]  [3] 

3.19-3 


The  problem  statement  implies  that  each  of  the  stubs  is  short-circuited  at  the  end  opposite 
the  main  line.  For  this  to  be  a  bandpass  filter,  the  magnitude  of  the  input  impedance  looking 
into  each  stub  must  be  very  high  -  nominally  infinite  -  since  then  the  filter  structure  would 
be  in  effect  the  main  line  by  itself,  with  no  stubs,  and  would  therefore  be  well-matched  at  the 
filter  input  and  output.  At  any  higher  or  lower  frequency  the  magnitude  of  the  stubs’  input 
impedance  can  only  be  less;  therefore,  the  input  impedance  of  the  filter  would  be  increasingly 
mismatched.  This  results  in  bandpass  response. 

The  shortest  length  for  which  the  magnitude  of  the  input  impedance  of  a  short-circuited 
transmission  line  is  infinite  is  A/4.  Therefore,  A/4  =  3.38  mm  and  subsequently  A  = 
13.52  mm  in  the  stub.  Therefore,  the  center  frequency  is 

f  _vp  _  0.6c 

;  “  X  “  ~T~ 


13.3  GHz 
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[m0091]  [4] 

3.19-4 


(a)  At  2.4  GHz,  the  free  space  wavelength  Ao  =  c//  =  12.5  cm.  Therefore  the  wave¬ 
length  in  the  line  A  =  0.67Ao  =  8.375  cm,  and  subsequently  the  length  of  each  section  is 
A/4  =  2.094  cm.  The  impedance  looking  into  each  stub  is  nominally  infinite  at  2.4  GHz; 
therefore  the  stubs  should  be  terminated  into  short  circuit  loads.  Then,  each  stub  will  trans¬ 
form  its  “load  impedance”  of  0  into  an  input  impedance  of  1/0  —>  oo  at  the  frequency  at 
which  it  is  a  quarter-wavelength  long.  The  resulting  design  is  shown  below: 


(b)  First,  note  that  the  input  impedance  of  a  short-circuited  stub  is  Zs  =  +jZotan/3Z  where 
Zq  is  the  characteristic  impedance  (50  H  in  this  case),  l  is  the  physical  length  of  the  stub 
(2.094  cm  in  this  case),  and  f3  =  2ir/\  (  =  27r//0.67c  in  this  case).  Consulting  Figure HTTl  we 
determine  that  the  response  at  a  specified  frequency  /  may  be  calculated  using  the  following 
steps: 


P 

Z, 

Z2 


Pl 

Pir. 


4~ 

4~ 

4~ 

4~ 

«- 


2  tt/ 

0.67c 
Zl  ||  Zs 
1  -I-  Tp-PP 

z°l  -  Te-iW  where  r  4 

z2 1|  Zs 

1  —  r  where  F  is  now 


Z\  —  Zq 
Z\  +  Zq 

A  Zin  Z  o 

Zin  +  Zq 


(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 


In  the  last  step,  Pl / Pm  is  response  as  defined  in  the  problem  statement.  This  expression 
works  under  the  assumption  of  no  loss  within  the  filter;  i.e.,  all  power  delivered  to  the  input 
is  subsequently  delivered  to  the  load,  and  none  is  dissipated  by  the  filter. 
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Frequency  Response  [Mag,  dB] 


Figure  3.1:  Schematic  representation  of  the  filter  with  a  matched  output  termination]. 


1  r  i  i  [  i  i 


-6  L 


1  1.5  2  2.5  3  3.5 

Freq  [GHz] 


42 


(c)  MATLAB  script  follows. 

clear  all; 
close  all; 

ZL  =  50.0;  7.  [ohm]  impedance  attached  to  output 
Z0  =  50.0;  °/„  [ohm]  characteristic  impedance 
1  =  0.02094;  %  [m]  section  length 

c  =  3 . Oe+8 ;  %  [m/ s] 

f.list  =  [1:0.01:3.8]* (le+9) ;  %  [Hz] 

n=0;  7o  counting  points 
for  f=f_list, 
n=n+l ; 

beta  =  2 . 0*pi*f / (0 . 67*c) ; 

Zs  =  +j*Z0*tan(beta*l) ; 

Z1  =  (ZL*Zs)/ (ZL+Zs) ; 

Gamma  =  (Z1-Z0) / (Z1+Z0) ; 

Z2  =  Z0*(l+Gamma*exp(-j*2*beta*l))/(l-Gamma*exp(-j*2*beta*l)) ; 
Zin  =  (Z2*Zs)/(Z2+Zs) ; 

Gamma  =  (Zin-Z0)/(Zin+Z0) ; 

P(n)  =  1-abs (Gamma) "2 ; 

end  l  for  f 

plot (f _list/(l . Oe+9) , 10 . 0*logl0 (P) , ;b-J ) ; 
grid  on; 

xlabelC’Freq  [GHz];); 

ylabel( ’Frequency  Response  [Mag,  dB] ’ ) ; 
axis ( [1  3.8  -6  +1]); 
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[m0090]  [1] 

3.20-1 


Summarizing  the  problem  statement:  P+v  =  5  W  and  Pl  =  4.6  W.  Therefore,  Pl/P^v  =  0-92. 
From  this,  we  may  find  the  magnitude  of  the  reflection  coefficient,  |r|,  using 


We  find  T  =  0.283  and 


SWR 


i  +  |r 
i  -  |r 


1.79 
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[m0090]  [2] 

3.20-2 


From  the  problem  statement,  F  =  0.3  +  j0.4  and  P+v  =  3  W.  Therefore, 


Pl 


(i  -  irf)  p+ 


2.25  W 
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[m0094]  [1] 

3.23-1 


From  the  problem  statement,  we  have  /  =  220  MHz,  antenna  impedance  Za  =  73  +  j'42 
and  characteristic  impedance  Z$  =  50  fl  for  both  the  transmission  line  and  the  stub.  The 
input  impedance  looking  into  a  length  d\  of  transmission  line  terminated  in  impedance  Za 

is 

1  +  Te~^dl 

Z1(W  =  Z.1_re_JW,  (3.18) 

where 

r  =  - -°  =  0.2719  +  jO.2486  (3.19) 

Z‘A  +  Z0 

The  first  task  is  to  find  the  smallest  /3d\  such  that  the  real  part  of  Yi (ftd\)  =  Z^ifddi)  equals 
Y0  =  Zq1  =  0.02  O-1.  After  a  few  minutes  of  trial  and  error  one  finds: 

=  1.345  rad)  =  0.0200  +  j'0.0159  fT1  (3.20) 

(You  could  also  do  this  with  a  Smith  chart  if  you  are  so  inclined.)  The  match  is  accomplished 
by  attaching  a  stub  having  input  admittance  Y2  =  — j0.0159  f)_1  in  parallel  with  Y\ ,  since 
then  the  combined  admittance  will  be  Y\  +  Y2  =  Y0  =  Zq1-  For  a  short-circuited  stub  of 
length  d2  we  would  want: 

Y2  =  —jY0  cot  j3d2  =  — j0.0159  f2-1  f3d2  =  0.900  rad  (3.21) 

For  an  open-circuited  stub  of  length  d2  we  would  want: 

Y2  =  +jY0ta,n/3d2  =  -j0.0159  H_1  =»  f3d2  =  2.471  rad  (3.22) 


The  short-circuited  stub  is  shorter,  so  that’s  the  preferred  solution.  All  that  remains  is  to 
figure  out  the  physical  lengths  from  the  electrical  lengths.  For  this,  we  need  to  know  (3.  The 
phase  velocity  is  vp  =  0.67c,  so 


2vr 


A  0.67  c/f 


=  6.8771  rad/m 


(3.23) 


Finally  we  have  the  solution: 


di  = 


/3d\  1.345  rad 

/3  6.877  rad/m 

/ 3d2  0.900  rad 


d2  — 


19.6  cm  distance  from  antenna  terminals  to  stub  (3.24) 

(3.25) 


/?  6.877  rad/m 


13.1  cm  stub  length 


and  the  stub  is  short-circuited 
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[m0094]  [2] 

3.23-2 


In  terms  of  the  variables  used  in  the  book,  the  problems  statement  is  indicating  that 
Y\  =  0.0128  —  jO. 0040  h2_1  and  that  Zin  is  real- valued.  Therefore,  Yin  =  1  / Zin  is  real- valued, 
and  must  be  equal  to  the  real  part  of  i.e. ,  Yin  =  0.0128  Q_1.  Therefore,  Zin  =  78.1  hi  , 
which  is  the  answer  to  part  (a). 


The  stub  is 
subsequently 


reing  used  to  cancel  the  imaginary  part  of  Yi,  so  Ystub 


Z stub  =  — j250  hi  ,  which  is  the  answer  to  part  (b). 


+j 0.0040  1  and 
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[m0094]  [3] 

3.23-3 


In  terms  of  the  variables  used  in  the  book,  the  problem  statement  indicates  Zl  =  35  —  j  10  hi, 
Zin  =  50  Q,  and  Z$  =  100  throughout.  The  voltage  reflection  coefficient  at  the  interface 
between  Zl  and  the  primary  line  is 


r 


Zl  ~  Z0 
Zl  +  Z0 


^  -0.473  -jO.  109 


Let  Y\  be  the  admittance  looking  into  the  primary  line: 


i  _  Yp-fiPh 

V  =  V 

1  °i  +  re-^i 


(3.26) 


(3.27) 


where 

Y0  =  =  0.01  mho  (3.28) 

Zq 

and  /3l\  is  the  electrical  length  of  the  primary  line.  To  match  the  real  part  of  the  admittances, 
we  require  Re  {17}  =  Re  {Ym}  where 


Yin  =  —  =  0.02  mho 

Z-J  in 


(3.29) 


is  the  input  admittance  corresponding  to  Zin.  Therefore  the  desired  value  of  (5l\  is  the 
solution  to 


Re 


1  -  Te~^h 
1  +  Te-^i 


2 


(3.30) 


Using  a  numerical  trial-and-error  search,  one  finds  j3li  =  0.362  rad.  Now  using  Equationl3.36l 


17  (fill  =  0.362  rad)  ^  0.0200  -  j0.0121  mho  (3.31) 

The  necessary  shunt  susceptance  (i.e. ,  the  imaginary  part  of  admittance)  is  —  Im  {17}  = 
+0.0121  mho,  since  this  will  cancel  the  susceptance  of  the  primary  line  when  placed  in 
parallel  with  the  primary  line.  Now  we  seek  the  shortest  stub  that  has  this  susceptance.  For 
an  open-circuited  stub  we  would  need 


+  lotan/3/2  =  +0.0121  mho  (3.32) 

where  /%  is  the  electrical  length  of  the  stub.  This  yields  =  0.8814  rad.  For  a  short- 
circuited  stub  we  would  need 


—  locot  pi2  =  +0.0121  mho  (3.33) 

This  yields  /3l2  —  —0.6893  rad.  /3  is  positive  and  length  can’t  be  negative,  so  we  need  the 
next  greater  value  of  /3l2  that  solves  the  above  equation.  Since  cot(-)  has  period  tt  radians, 
the  desired  value  is  2.4522  rad.  This  is  much  longer  than  the  result  for  the  open-circuited 
stub,  so  we  choose  the  open  circuit  result. 
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(a)  Note 

pi  =  21. 1  =  2tt 
A  A 

Therefore  to  express  electrical  length  in  wavelengths,  we  simply  divide  by  2n.  Thus,  the 
solution  to  the  problem  is: 

Primary  line  length  l\  =  0.058A 
Stub  length  l2  =  0.140A 
Stub  is  open-circuited. 

(b)  Note 

P  =  —  =  ^  48.332  rad/m 

vp  0.65c 

since  /  =  1.5  GHz  and  the  velocity  factor  is  65%.  Therefore 

,  Ph  ^  - 

<i  =  — —  =  7.5  mm 

1 9  =  =  18.2  mm 

P  - 


(c)  Let’s  define  T*n  as  the  voltage  reflection  coefficient  at  the  input  of  the  matching  structure. 
(Note  that  this  is  different  from  T  defined  in  previous  parts,  which  is  the  voltage  reflection 
coefficient  at  the  output  of  the  matching  structure.)  Therefore  the  fraction  of  power  delivered 
(Pl)  to  power  incident  (P+,)  is: 


Pl 

PaX 


(3.34) 


where 


Zin  -  50  Q 

zin  +  so  n 


(3.35) 


and  where 

/  u  Te-j2f}l1  \ 

Zin  =  (-jZ0  cot  pip  ||  ^Z0  -  _  Te_mi  J  (3.36) 

Be  careful:  The  sweep  in  frequency  appears  as  a  sweep  in  the  value  of  P  in  the  above  equa¬ 
tion.  A  plot  of  the  result  follows. 
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Power  Delivered  /  Power  Incident  [linear  units] 


frequency  [GHz] 
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Chapter  5 
Electrostatics 
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[m0102]  [1] 

5.1-1 


The  electric  field  due  to  a  point  charge  q  is 


E(R) 


R 


g 

A'KeR2 


where  R  is  the  position-free  vector  pointing  from  the  charge  to  the  field  point.  From  the 
problem  statement,  q  =  —24  nC,  er  =  2,  and 


R  =  x  +  y2  +  z3  m 


Thus 


Thus 


R  =  |R|  =  \/l2  +  22  +  32  =  3.74  m 

R. 

R  =  —  =  x0.267  +  y0.534  +  z0.802 

it 

e  =  ere0  =  2  •  8.854  x  10-12  F/m 


E(R) 


— x2.06  —  y4.12  —  z6.18  V/m 
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[m0103]  [1] 

5.2-1 


From  the  problem  statement, 

q\  =  +3  nC  at  ri  =  — z d  and 

g2  =  +3  nC  at  r2  =  +z d  where  d  =  0.5  m; 

the  held  point  of  interest  is  r  =  +xx  where  x  =  +1.5  m,  and 

e  =  e0  =  8.854  x  10~12  F/m. 


The  electric  held  intensity  due  to  two  point  charges  is 


E(r)  = 


47T6 


ri 


i9i 


r2 


r?2 


ri 


r2 


In  this  problem: 


r  —  ri  =  +xx  —  (— z  d)  =  +xx  +  z  d 


ri 


=  Vx2  +  d2 


r  —  r2  =  +xx  —  (+z  d)  =  +xr  —  z  d 


f2 


=  V  x2  +  d2 


Substituting: 


E  (+xx)  = 


47ren 


+xr  +  z  d  +xa;  —  z  d 

r9+  — : - 


(x2  +  d2f2  (: x 2  +  d2f2 


where  we  have  made  the  definition  q  —  qi  —  qi-  Note  that  the  z- directed  components 
cancel,  as  expected  from  the  symmetry  of  the  problem.  Eliminating  these  components  and 
simplifying: 


E(+xx)  =  x- 


q 


x 


2vre0  (a;2  +  rf2)3/2 

Now  take  a  moment  to  conhrm  that  the  solution  is  dimensionally-correct  and  makes  physical 
sense.  Finally,  substituting  values,  we  obtain: 


E(+xl.5  m)  =  +x(20.5V/m) 


For  a  single  charge  go  at  the  origin  to  create  this  held,  we  require 


qo 


47re0xs 


=  +x  (20.5  V/m) 


which  yields  go 


+5.12  nC 
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[mOlOO]  [1] 

5.3-1 


From  the  problem  statement,  pv  =  Kr~2  where  K  =  2  C/m.  From  dimensional  analysis, 
it  is  clear  that  this  is  a  volume  charge  density.  We  seek  the  total  charge  Q  in  a  volumetric 
region  V  bounded  by  the  constant-coordinate  surfaces  r  =  a  and  r  =  b  where  a  —  1  m  and 
b  =  2  m.  In  general, 

Q  —  /  pvdv 
Jv 

In  this  case,  using  spherical  coordinates: 


Q  = 


r'2n 


’  r=a  J  6=0  J  6=0 


(r2  sin  6  dr  d6  d(j)  ) 


Factoring  into  separate  integrals: 


-  rb 

/*7T 

r 

Q  =  K 

/  dr 

/  sin  6  d6 

/  d(f> 

_J  r=a 

Jo= 0 

J<p= 0 

Evaluating  the  integrals: 


Q  =  K  [b  -  a]  [2]  [2tt]  =  4ttA'  ( b  -  a) 

This  a  good  point  at  which  to  check  for  dimensional  consistency  (i.e.,  correct  units). 


Using  the  given  values  of  K,  a,  and  b  we  obtain  Q  = 


25.1  C 
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[mOlOO]  [2] 

5.3-2 


From  the  problem  statement,  the  volume  charge  density  is 

_  Pv  0 

Pv  o  •  r\ 

rz  sm  9 


(5.1) 


for  r\  <  r  <  r2,  where  r±  =  1  m  and  =  2  m;  and  6\  <  6  <  62,  where  6\  =  7t/4  (rad)  and 
62  =  37t/4.  Also  pv 0  =  1.3  C/111  inside  these  limits,  and  pv 0  =  0  outside  these  limits.  Let  V 
be  the  region  of  space  where  pv0  7^  0.  Then  the  total  charge  Q  is 


=  Pv 0  {r 2  -  t  1)  (6*2  -  9 1)  (27t) 


(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 


Note  this  result  is  dimensionally  correct.  Substituting  the 


tain  Q  =  12.83  C 


values  established  above,  we  ob- 
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[m0104]  [1] 

5.4-1 


Interpreting  the  problem  statement: 

pSi i  =  +4  nC/m2  for  sheet  in  the  x  =  0  plane.  Let  the  field  from  this  sheet  be  Ei. 
pS) 2  =  +16  nC/m2  for  sheet  in  the  y  =  0  plane.  Let  the  field  from  this  sheet  be  E2. 
pSj 3  =  +64  nC/m2  for  sheet  in  the  z  —  0  plane.  Let  the  field  from  this  sheet  be  E3. 
Also,  e  =  2e0. 


The  electric  field  intensity  due  to  a  single  sheet  of  charge  having  charge  density  ps  in  the 
z  —  0  plane  is  worked  out  in  the  book.  It  is: 

-  Ps 

+z— sgn  2 

This  corresponds  to  the  third  sheet  of  charge  above.  Since  the  region  of  interest  is  z  >  0: 


Similarly, 

E,  = 

4e0 

E2  =  +y^ 
4e0 

The  total  field  is  the  sum  of  these  three  fields.  Thus: 


-tp  -  Ps,  1  .  *  Ps, 2  .  Ps, 3 

E  =  x- - h  y-: - h  z— J— 

4e0  4e0  4e0 


Substituting  values,  we  obtain: 

E  = 


x  (113  V/m)  +  y  (452  V/m)  +  z  (1807  V/m) 
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[m0104]  [2] 

5.4-2 


From  the  problem  statement,  we  have  electric  field  intensity  E/me  due  to  line  charge  density 
pi  =  +8  mC/m  along  the  2-axis.  Thus, 

E  iine  =  P^~  (5.7) 

ZTTep 

Also,  we  have  electric  field  intensity  E Sheet  due  to  surface  charge  density  ps  =  +12  mC/m2 
in  the  2  =  0  plane.  Thus, 

E sheet  =  for  2  >  0  (5.8) 

The  total  electric  field  is  determined  by  superposition: 

E  =  EUne  +  Esheet  =  p-^—  +  for  2  >  0  (5.9) 

27rep  2e 

Note  that  this  is  dimensionally  correct.  Also  from  the  problem  statement  we  have  e  =  ereo 
where  eo  =  8.854  x  10-12  F/m  and  er  =  2.  Finally: 


E 


A  71.9  MV  A 

p - h  Z 

P 


^338.8 


for  2  >  0 


(5.10) 
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[m0014]  [1] 

5.5-1 


Note  F12  =  Q2E1  where  E|  is  the  electric  held  intensity  associated  with  Q1.  Thus: 

F  =  —  -  R  Q1 
1  Q2  1247rel?f2 

We  can  write  this  in  terms  of  the  electric  flux  density,  assuming  a  isotropic  and  homogenous 
medium: 

Now  let’s  put  Q 1  at  the  origin,  and  let  S  be  a  sphere  of  radius  a  centered  at  the  origin. 
Then,  the  left  hand  side  of  Gauss’  Law  is: 


/*7T  /»27T 

D  ■  ds  =  /  / 

s  Je= 0  J<j)= 0  . 


.  Q 1 
47ra2 


/*7r  r2n 

r™2sin0  dO  d<f>  1  =  —  /  /  sin  9  dO  d<f>  =  Q 

J  4?r  i<*=o 


Because  we  put  Q 1  at  the  origin  and  defined  5  to  surround  it,  Qenci  =  Qi,  which  is  what  we 
expect  from  Gauss’  Law.  Therefore,  Coulomb’s  Law  is  a  solution  to  -  a  special  case  really 
-  of  Gauss’  Law. 


You  might  be  inclined  to  object  on  the  grounds  that  Gauss’  Law  doesn’t  say  anything  about 
force  or  electric  held  intensity.  This  is  true!  However,  electric  held  intensity  is  defined  by 
force;  i.e. ,  E|  =  F12/Q2  is  a  definition  for  E1;  and  not  derived  from  something  else.  Similarly, 
Dx  =  eEi  is  a  definition  for  D1;  and  not  derived  from  something  else.  So,  Gauss’  Law  is  as 
fundamental  as  it  gets. 


[m0014]  [4] 

5.5-2 


By  symmetry,  there  can  be  no  variation  in  the  z  or  0  dimensions.  Therefore,  our  answer 
can  depend  only  on  p.  The  three  regions  to  consider  are  inside  the  inner  surface  of  the  shell 
(p  <  1  m),  inside  the  shell  itself  (1  <  p  <  3  m),  and  outside  the  outer  surface  of  the  shell 
(p  >  3  m). 

The  integral  form  of  Gauss’  Law  is: 


D  •  ds  =  Q 


encl 


(5.11) 


where  S  is  any  closed  surface.  Since  we  are  asked  for  electric  held  intensity,  we  can  use 
D  =  eE  to  obtain: 


E  *  ds  Q  end 

e 


(5.12) 


Note  that  we  use  e  as  opposed  to  e0,  since  the  latter  infers  free  space  conditions,  and  we 
haven’t  been  told  that. 


For  p  <  1  m,  Qenci  =  0;  i.e.,  there  is  no  surface  that  we  can  define  that  encloses  charge. 
Combined  with  the  symmetry  argument,  we  have  that 


E  =  0  in  this  region. 


For  1  <  p  <  3  m,  Qenci  depends  on  p.  Combined  with  the  symmetry  argument,  we  have 
E(p)  =  pE(p).  Thus,  a  good  choice  for  S  is  a  cylinder  centered  on  the  0  axis.  This  gives  us: 


r+L/2 


'  z=—L/2 


pE(p)  ■  p  p  d(j)  dz  +  0 


r*+L/2 


'  Z——L/2 


pv  dp  pd(j)  dz 


(5.13) 


where  “+0”  on  the  left  hand  side  is  the  contribution  from  the  constant-^  surfaces  (the  “end 
caps” )  of  the  cylinder  -  zero  because  the  normal  to  those  surfaces  (z)  is  perpendicular  to  E. 
Now  evaluating: 


E(p)  2'KpL 


-  pvnL(p2  -  1  ill2) 


(5.14) 


Finally: 


E(p)  =  p 


Pv 

2e 


in  this  region. 


(5.15) 


Assuming  pv  is  in  C/m3  and  e  is  in  F/m,  E(p)  will  be  in  V/m.  However,  if  you  say  simply 
“1”  as  opposed  to  “1  m2”  in  the  above  expression,  then  you  must  indicate  the  units  of  p 
(being  meters)  as  well.  In  electromagnetics,  a  powerful  (but  unappreciated)  technique  for 
checking  your  work  is  to  make  sure  your  solution  has  the  right  units.  This  is  called  dimen¬ 
sional  analysis.  You  should  be  able  to  substitute  units  for  each  of  the  quantities  in  the  above 
solution  and  find  that  the  result  has  units  of  V/ m  -  can  you  do  this? 
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For  p  >  3  m,  Qenci  is  constant  at  the  maximum  value  (since  all  the  charge  has  been  enclosed), 
but  E(p)  is  still  a  variable  function  of  p.  Thus,  the  left  hand  side  of  Equation  15. 141  remains 
the  same,  but  the  right  hand  side  is  evaluated  at  p  =  3  m.  This  yields: 


E(p)  2-KpL 


-  PvTtL  ■  8  m2  . 
e 


(5.16) 


Thus: 


E(p)  =  p 


Pv  4  m2 
e  P 


in  this 


region. 


(5.17) 


Suggestion:  Try  dimensional  analysis  on  this  solution.  Can  you  see  why  it  is  important  to 
say  “4  m2”  as  opposed  to  just  “4”? 
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[m0014]  [2] 

5.5-3 


By  symmetry,  there  can  be  no  variation  in  the  9  or  0  dimensions.  Therefore,  our  answer 
can  depend  only  on  r.  The  three  regions  to  consider  are  inside  the  inner  surface  of  the  shell 
(r  <  2  m),  inside  the  shell  itself  (2  <  r  <  4  m),  and  outside  the  outer  surface  of  the  shell 
(r  >  4  m). 

The  integral  form  of  Gauss’  Law  is: 


D  •  ds  =  Q 


enclosed 


(5.18) 


where  S  is  any  closed  surface.  Since  we  are  asked  for  electric  held  intensity,  we  can  use 
D  =  eE  to  obtain: 


E  •  ds  Q enclosed 

e 


(5.19) 


Note  that  we  use  e  as  opposed  to  e0,  since  the  latter  infers  free  space  conditions,  and  we 
haven’t  been  told  that. 


For  r  <  2 
Therefore, 


m,  Q enciOSed  =  0;  i.e.,  there  is  no  surface  that  we  can  define  that  encloses  charge. 


E  =  0  in  this  region. 


For  2  <  p  <  4  m,  Qenciosed  depends  on  r.  Combined  with  the  symmetry  argument,  we  have 
E(r)  =  r E(r).  Thus,  a  good  choice  for  S  is  a  sphere  centered  at  the  origin.  This  gives  us: 


r*27T 


r E(r)  ■  r  r2  sin  9  dO  d<f> 


t>= o 


-2  7T 


C  .1  T  —  O  ./  ff  —  l 


pv  r 2  sin  9  dr  d9  dcj) 


r= 2  J  6= 0  Jd>= 0 


(5.20) 


Now  evaluating: 


Finally: 


Airr2  E(r)  =  —  - —  (r3  —  [8  m3]) 
6  3 


in  this  region. 


(5.21) 


(5.22) 


Assuming  pv  is  in  units  of  C/m3  and  e  is  in  F/m,  E(p)  will  be  in  V/m.  However,  if  you 
say  simply  “8”  as  opposed  to  “8  m3”  in  the  above  expression,  then  you  must  indicate  the 
units  of  r  (being  meters)  as  well!  In  electromagnetics,  a  powerful  (but  sadly,  unappreci¬ 
ated)  technique  for  checking  your  work  is  to  make  sure  your  solution  has  the  right  units. 
This  is  called  dimensional  analysis.  You  should  be  able  to  substitute  units  for  each  of  the 
quantities  in  the  above  solution  and  find  that  the  result  has  units  of  V/m  -  can  you  do  this? 


Also  note  that  you  have  a  second  way  to  check  your  solution  -  it  must  be  equal  to  the 
solution  for  the  first  region  for  r  =  2  m.  Note  that  it  is. 
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For  r  >  4  m,  Qenciosed  is  constant  at  the  maximum  value  (since  all  the  charge  has  been 
enclosed),  but  E(r)  is  still  a  variable  function  of  r.  Thus,  the  left  hand  side  of  Equation  15.211 
remains  the  same,  but  the  right  hand  side  is  evaluated  at  r  =  4  m.  This  yields: 


47rr2  E(r)  —  —  —  (56  m3)  . 
e  3 


(5.23) 


Thus: 


E(r) 


„  pv  56  m3 
r  3e  r2 


in  this  region. 


(5.24) 


Suggestion:  Try  dimensional  analysis  on  this  solution.  Can  you  see  why  it  is  important 
to  say  “56  nr3”  as  opposed  to  just  “56”?  Also,  confirm  that  your  answer  agrees  with  the 
Region  2  answer  for  r  =  4  m. 
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[m0014]  [3] 

5.5-4 


a)  A,  B,  and  C  have  units  of  V/m4,  V/m3,  and  V/m2  ,  respectively. 


(b)  According  to  the  integral  form  of  Gauss’  Law: 


Qend  =  f  D  •  ds  =  e0  f  E  •  ds 
Js  Js 

Here,  S  is  the  surface  of  the  box-shaped  region,  and  ds  is  the  normal  to  each  of  the  six  sides. 
This  integral  is  easiest  to  handle  as  the  sum  of  integrals  over  each  side,  since  then  ds  will  be 
constant  over  each  of  these  integrals.  Here  we  go: 

/ry=+ 1  rz= o  2 

E  •  (—x  dy  dz)  —  —  I  /  A  (—1)  z 2  dy  dz  =  +-d 
-x  side  Jy= o  Jz=- i  3 

**/=+ 1  rz= 0 


r  ry=+i  /•*= u 

/  E  ■  (+x  dy  dz)  —  I  A  (+1)  z2  dy  dz  =  +-d 

J+x  side  Jy= o  Jz=- i  3 

/px=-\- 1  rz= 0 

E  •  (— y  dx  dz)  =  —  /  (—5  (0)  z)  dx  dz  —  0 

-y  side  Jx=-i  Jz=-i 

r  fX=  + 1  nz= 0 

/  E  •  (+y  dx  dz)  —  I  /  (—5  (+1)  2)  dx  dz  =  +B 

'+y  side  Jx=-\  Jz=- 1 

E  •  (— z  dx  dy)  =  —  /  C'x  dx  dy  =  0 

-z  side  Jx=- 1  Jy=o 

r  rx=+l  py=- 1-1 

/  E  ■  f+z  dx  dy)  =  /  Cx  dx  dy  =  0 

J+z  side  Jx=-i  Jy= 0 


So  we  find: 


Qencl  —  e0  (  0^4  +  -A  +  0  +H  +  0  +  0  1—  Co  (  —A  + 


having  units  of  Coulombs  if  e0  is  in  F/m  and  the  dimensions  are  all  in  meters.  The  chances  for 
units-related  confusion  is  reduced  if  consider  what  has  actually  happened  in  the  integration 
and  say  specifically: 


Q 


encl 


eo 


A  +  (1  m4)  B 
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[m0014]  [5] 

5.5-5 


The  problem  is  easily  solved  using  Gauss’  law  in  integral  form: 

Qenci  —  j)  D  •  ds  (5.25) 

where  S  is  any  surface  which  completely  surrounds  the  charge,  D  is  the  electric  flux  density, 
and  ds  is  the  differential  surface  element.  The  easiest  surface  in  this  case  is  a  sphere  of  radius 
tq,  centered  on  the  origin,  with 


r0  > 


Vl2  +  l2  +  l2 

2 
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2 


(5.26) 


Note  that  it  is  not  important  for  the  radius  of  the  sphere  to  be  close  to  this  number;  it  is 
merely  necessary  that  the  radius  be  greater  than  this  number.  In  fact,  we  shall  see  below 
that  the  radius  doesn’t  matter  at  all,  as  long  as  it  is  at  least  this  big. 


From  the  problem  statement: 


E  =  r 


3  V  •  m 


(5.27) 


The  problem  indicates  the  medium  is  free  space,  so  the  permittivity  e  =  eo  =  8.854  x 
10-12  F/m.  Therefore: 

_  „  3  V  •  m 

D  =  e0E  =  re0 5 — 


Now  putting  this  all  together: 


Qenci 


r>7r  n2n 


'  9=0  J(f>= 0 


re0- 


3  V  •  111 


)■< 


r  Tq  sin  6  d6  d<fdj 


9=0  J<t>= 0 


=  e0  (3  V  •  m)  f 
Je. 

=  e0  (3  V  •  m)  f  f  sin  9 

\Je= 0 

=  e0  (3  V  •  m)  (2)  (27t) 


sin  0  d6  d(j) 


d9 


r*27T 


d(p 


1 4>= 0 


(5.28) 

(5.29) 

(5.30) 

(5.31) 

(5.32) 


Note  this  result  is  dimensionally  correct.  Substituting  the  values  established  above,  we  ob¬ 


tain  Qend  =  333.8  pC 
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[m0149]  [1] 

5.6-1 


This  is  essentially  the  same  problem  shown  as  an  example  in  the  book,  for  which  the  electric 
held  intensity  was  found  to  be 


E  =  p 


Pi 

27rep 


where  here  pi  =  —2.1  mC/m  and  p  is  the  distance  from  the  z-axis.  The  electric  flux  density 
is  D  =  eE,  so  the  permittivity  doesn’t  matter.  The  result  is: 


D  =  p 


Pi 

2ttp 


-P 


334  /iC /m 

P 
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[m0045]  [1] 

5.7-1 


From  Gauss’  Law,  pv  =  V  •  D  =  V  ■  (e0erE)  =  e0erV  ■  E.  Calculating  the  divergence: 


v ' E  =  £  v/m2^) +  §y  (( 2  v/m3M  +  §z  ((x  v/m 3M 

=  (6  V/m2)(l)  +  (2  V/m3)(z)  +  (1  V/m3)(0) 

=  (6  V/m2)  +  (2  V/m3)z 
e0  =  8.854  x  1CT12  F/m  and  er  =  4.5,  so 


Pv 


239.1  pC/m3  +  (79.7  pC/m 4)z 
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[m0045]  [2] 

5.7-2 


(a)  pv  =  V  ■  D  =  e0V  •  E.  Here  it’s  easiest  to  use  Cartesian  coordinates,  for  which 

_  A  d  „  d  A  d 

V  =  +yIT  +  zb~~ 
ox  ay  oz 


so  we  have 


Pv  £0 


d 


d 


d 


*dx+ydy+Zdz 


[x  (2  V/ m)  sin  x  cos  y  —  y  (2  V/ m)  cos  x  sin  y\ 


=  e0  [  (2  V/ m2)  cos  x  cos  y  —  (2  V/ m2)  cos  x  cos  y\  =  _0_ 

This  is  an  example  of  a  divergence- free  field.  It  seems  that  there  can  be  an  electric  field  even 
when  there  is  no  charge.  This  means  simply  that  the  source  charge  must  lie  entirely  outside 
the  region  begin  considered. 


(b)  In  this  case  we  have 


Pv 


d 


d 


d 


*dx  +ydy  +Zdz 


[x  (3  V/ m)  cos  xy  +  y  (3  V/ m)  sin  xy\ 


=  (3  V/ m3)  Co  [—y  sin  xy  +  x  cos  xy] 
which  has  the  expected  units  of  C/m3. 
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[m0045]  [3] 

5.7-3 


(a)  A,  B ,  and  C  have  units  of  V/m4,  V/m3,  and  V/m2  ,  respectively. 


(c)  According  to  the  differential  form  of  Gauss’  Law,  we  have  for  this  problem: 


pv  =  e0V  •  E  =  e0 


A(Ax^  +  |-( -Byz )  +  §-z(Cx) 


=  e0  [ Az 2  —  Bz  +  0]  . 


This  is  charge  density  as  a  function  of  position.  The  enclosed  charge  is  obtained  by  integrat¬ 
ing  over  the  region  of  interest: 


Q  encl  /  Pv  dv 


f'X=-\-l  ry=-\-l  nz= 0 


'V 


'  x=—l  J  y= 0  J  z=— 1 


e0  (Az2  —  Bz)  dx  dy  dz 


The  integrations  over  x  and  y  factor  out  and  are  equal  to  2  and  1  respectively.  What’s  left 
is: 


Qenci  —  (2)(l)eo 


-z=+ 1 


'  z=- 1 


(Az2  —  Bz)  dz  = 


Qend  CO 


-  nr 
3 


^  A  +  (1  m4)  B 


Note  that  this  result  agrees  with  the  result  obtained  using  the  more  direct  approach  of  us¬ 
ing  the  integral  form  of  Gauss’  Law.  You  should  note  that  the  reason  the  results  are  the 
same  is  not  really  related  to  electromagnetics,  but  rather  due  to  the  Divergence  Theorem 
(of  mathematics),  which  relates  the  behavior  of  a  vector  field  in  a  volume  to  the  behavior  of 
that  same  vector  field  over  the  enclosing  surface. 


[m0061]  [1] 

5.8-1 


The  change  in  the  energy  of  the  system  resulting  from  moving  the  particle  a  small  distance 
IAI  is: 

AW  «  -qE  ■  IAI 

Power  P  is  energy  per  time,  so  the  power  required  to  do  this  is: 


% 


AW 

At 


— qE  ■ 


IAI 
A t 


where  At  is  the  time  required  for  the  particle  to  traverse  the  distance.  Note  that  in  the 
limit  as  At  — >  0,  \Al/At  is  the  velocity  v  of  the  particle.  Taking  the  limit  and  making  the 
substitution, 

P  =  — gE  ■  v 

This  is  the  “instantaneous  power”  required  at  time  t  and,  through  t,  the  position  r (t). 


Interpreting  the  problem  statement:  q  =  —4  mC; 

E  =  E0 z,  where  E0  =  3  V/m;  and 

r(t)  =  xacostnZ  +  ybsinut  +  z ct,  where  a  =  b  =  2  m,  oj  =  n  rad/s,  and  c 
Note 

Therefore, 


d  „ 

— r (t)  =  —xacu  sm cot  +  ybu  coscot  +  zc 

ill 


4  m/s. 


P  =  —qE  ■  v  =  —qE0c 


(—4  mC)  (3  V/m)  (4  m/s) 


48  mW 
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[m0064]  [2] 

5.12-1 


Let  us  arbitrarily  assume  the  charge  is  aligned  along  the  z  axis.  Then  the  electric  field 
intensity  is  given  by 


E  0)  =  p 


Pi 


2-Kep 


(5.33) 


The  potential  difference  is: 


V21  =- 


/•point  2 
'point  1 


rp  2 


Pi 


E  •  d\  —  —  /  p- - •  pdp  = 


'pi 


2nep 


In  — 
27 re  p2 


(5.34) 
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[m0064]  [1] 

5.12-2 


The  electric  field  intensity  resulting  from  a  single  line  of  uniform  charge  density  is  given 
by 


E(R) 


R 


Pi 

2neR 


(5.35) 


where  R  is  the  distance  between  the  point  of  interest  and  closest  point  on  the  line,  R  points 
from  that  point  on  the  line  to  the  point  of  interest,  and  R  =  HR. 


If  we  have  two  such  lines  of  charge,  then  by  superposition  we  could  write: 


E(R) 


Ri 


Pi 

2'KeRi 


+  R-2 


Pi 

27Tei?2 


(5.36) 


where  the  subscripts  “1”  and  “2”  refer  to  the  geometry  relative  to  the  first  and  second  lines 
of  charge,  respectively.  In  this  problem,  the  lines  of  charge  and  the  point  of  interest  all  lie 
in  the  x  —  y  (z  —  0)  plane.  Since  this  is  the  case  we  may  write  simply: 


E(x,y)  =  x-^-  +  y-^~  for  £  =  0.  (5.37) 

27 rex  2-Key 

The  potential  difference  is: 

/•point  2 

V21  =  -  /  E  •  dl  .  (5.38) 

./point  1 

Remember  that  the  answer  should  be  the  same  for  any  path  between  the  points,  so  you 
might  as  well  choose  one  that  makes  the  problem  simple.  Here’s  the  result  using  one  of  two 
equally-easy  paths: 


E2i  = 


-  Pi  ,  *  Pi 

x - b  y - 

2irex  27re|/2 


Pi  ^  Pi 

2k£X\  J  27767/ 


•  ydy 


(5.39) 


That  is,  first  move  from  (2,4)  m  to  (1,4)  m  along  the  y  —  y2  —  4  m  line,  and  then  move 
from  (1,4)  m  to  (1, 1)  m  along  the  x  —  X\  —  1  m  line.  Evaluating: 


V21  = 


Pi 

2776 


0.331  - 
6 


(5.40) 


Here,  you  can  check  your  results  using  dimensional  analysis  (C/m  divided  by  F/m  gives  C/F 
=  V).  You  can  also  check  that  the  sign  is  correct:  Point  2  is  closer  to  both  lines  of  charge 
than  point  1,  so  when  the  charge  is  positive,  work  is  being  done  and  the  potential  difference 
is  positive.  Said  differently,  the  potential  at  point  2  is  higher  than  the  potential  at  point  1. 
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[m0064]  [3] 

5.12-3 


The  point  charge  q\  =  +3  /iC  at  the  origin  creates  a  potential  field 

V(r)  =  £  (5.41) 

where  r  is  the  distance  from  the  origin.  Since  the  permittivity  of  the  medium  is  specified  to 
be  twice  that  of  free  space,  e  =  ere0  =  2e0.  Thus: 


V(r) 


‘h 

87re0r 


(5.42) 


The  potential  difference  V2i  at  r2  relative  to  r!  is  independent  of  the  path  taken  between  the 
points;  it  depends  only  on  the  endpoints.  Thus: 


where  r  1 


v21  =  v (r2)  -  y(n)  = 

8vre0r2 

^/32  +  (— 4)2  =  5  m  and  r2  =  1  m. 


gi  gi  fl 

8vre0r  1  87re0  \r2  ’’I/ 


Therefore  V21 


+10.8  kV 


(5.43) 


Note  that  the  result  does  not  depend  on  the  value  of  the  charge  (g2  =  +2  /xC)  being  moved 
from  ri  to  r2.  This  is  the  whole  point  in  defining  a  scalar  electric  potential:  It  describes 
energy  in  the  field  independently  of  the  charge  that  experiences  it.  If  necessary,  one  may 
subsequently  calculate  the  energy  associated  with  this  potential  difference  as  q2V i2. 
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[m0063]  [1] 

5.14-1 


To  begin,  it  will  be  convenient  to  first  convert  r0  from  Cartesian  to  spherical  coordinates. 
Here  we  go: 


r o  =  \Jxl  +  yl  +  zl  =  5.38  cm 
0O  =  arccos  ( z0/r0 )  =  42.0° 
0o  =  arctan  (y0/xo)  =  56.3° 


V (r0)  =  Vorp  cos  6 o  =  10.8  mV 


(b) 


E  =  -VH  =  -[r|-  +  0-4  +  O]H 
or  r  89 


The  third  term  in  the  gradient  is  zero  because  V  in  this  problem  does  not  vary  with  0. 
Continuing: 

^dV  ~ldV 

E  =  — r— - 9-—-  =  —r2V0r  cos 9  +  9V0r  sin 9 

or  r  89 

So:  _ 


E(r0)  =  — r400  +  0180  mV/m 


(c) 


Pv  —  V  •  D  —  6qV  •  E  —  6q 


^4  (' r2Er )  +  ~ TTlIi  (Eo  sin0)  +  0 


r  sin  9  89 


where  Er  and  Eg  are  the  r-  and  0-directed  components  of  E.  The  third  term  in  the  divergence 
is  zero  because  E  in  this  problem  does  not  vary  with  0.  Note: 

1  c)  1  d  1  ^  1 

(r2Er)  =  —  —  (r2  (—2V0r  cos  9))  =  (~2V0r3  cos  9)  =  —  (-6Vor2cos0)  =  -6Vocos0 


ry  2 


2 

and 

1  ^  |  ^  |  ^  | 

— —77  (Eg  sin  0)  =  — 7  —  "7T7T  ((Vor  sin0)  sin0)  =  — -77  (Vorsin20)  =  — —  (2 V0r  sin  0  cos  0) 
r  sm  0  <70  r  sm  0  a0  r  sin  0  o0  r  sm  0 


=  2Vq  cos  0 


Continuing: 


Pv  =  Co  [ — 6V0  cos  0  +  2Vo  cos  0]  =  — 4Voe0  cos  0 
(Good  time  for  a  units  check...)  At  the  point  of  interest: 


Pv(tq)  —  -131  pC/1113 
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[m0063]  [2] 

5.14-2 


From  the  problem  statement, 
where  Vq  =  4  V  ■  m1/2.  So: 

d 


—  TZ.-r-V2 


V(r)  =  V0r 


E{  r)  =  —  W  (r) 


=  — r— Vo r  +  terms  that  go  to  zero  because  ^ 


dr 


d_ 

d<j> 


=  — rVh  (  — r  3/2 


Vr 


=  +r — r 
2 


°  -3/2 


+r  (2  V  •  m1/2)  r_3/2 


Note  that  the  answer  is  dimensionally  correct  (and  unambiguously  so). 


0 
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[m0067]  [1] 

5.15-1 


Poisson’s  Equation  is 

V2E  =  -P—  .  (5.44) 

e 

The  geometry  of  the  problem  suggests  cartesian  coordinates,  and  symmetry  such  that  dV/dy  = 
dV/dz  =  0  is  implied.  Thus,  the  above  equation  becomes: 


dx 2 


V(x) 


Pvix) 

e 


Integrating  both  sides  with  respect  to  x  we  have 


d_ 

dx 


V(x) 


pv(x)  dx  . 


(5.45) 


(5.46) 


Integrating  both  sides  again  with  respect  to  x  we  have 


V(x) 


pv(x)  dx 


— OO  LJ  —  OO 


dx  . 


(5.47) 


The  first  chore  is  to  take  care  of  that  sequence  of  integrations.  We  begin  with  mathematical 
restatement  of  the  given  volume  charge  density: 


{0  ,  —  oo  <  x  <  —  b 

—a  ,  —  b  <  x  <  0 
+a  ,  0  <  x  <  +b 

0  ,  +  b  <  x  <  —oo 


Integrating  once: 


pv(x)dx 


0  ,  —  oo  <  x  <  —  b 

—a(x  +  5),  —  b  <  x  <  0 
+ax  —  ab  ,  0  <  x  <  +b 

0  ,  +  b  <  x  <  —oo 


(5.48) 


(5.49) 


If  you  have  a  hard  time  seeing  this,  consider  sketching  pv(x)  and  then  doing  the  integration 
graphically.  Integrating  the  second  time: 


pv(x)dx 


—  OO  L*'  —  OO 


dx 


0 

—(a/2)x2 

+{ci/2)x2 

-ab2 


abx  —  ab2/ 2  , 
abx  —  ab2  / 2  , 


—  oo  <  x  <  —b 
—  b  <  x  <  0 
0  <  x  <  +b 
+  b  <  x  <  —oo 


Substituting  this  into  Equation  15.471  we  obtain: 


{0  ,  —  oo  <  x  <  —  b 

+  (a/2)x2  +  abx  +  ab2/ 2  ,  —  b  <  x  <  0 
—  (a/2)x2  +  abx  +  ab2/ 2  ,  0  <  x  <  +b 

ab2  ,  +  b  <  x  <  —oo 


(5.50) 


(5.51) 
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Now  check  your  answer.  First,  as  always,  check  that  is  is  dimensionally  correct.  Second, 
note  that  V(x)  should  be  a  continuous  function  of  x ,  since  integration  over  any  function 
(specifically  excluding  the  impulse  or  “delta”  function)  results  in  a  continuous  function. 


To  find  the  volume  charge  density  a  in  terms  of  V,  we  simply  evaluate: 


Vd  =  V(x  =  +b )  -  V(x  =  - b ) 


and  solve  for  a: 


e 


ab 2 
e 


(5.52) 


(5.53) 


Given  the  relative  permittivity  of  silicon  er  fa  12,  b  —  100  // m ,  and  Vd  =  0.4  V,  we  find 
a  =  4.25  mC/m3  . 
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[m0067]  [2] 

5.15-2 


The  symmetry  of  this  problem  suggests  a  solution  in  spherical  coordinates.  Laplace’s  Equa¬ 
tion  in  spherical  coordinates  is: 


1  d 


V2V  =  —  —  r2——  + 


,dV\ 


d 


r2  dr  \  dr  J  r2  sin  0  dO 


n  dV\ 

sill  6»  — —  + 


d2V 


d6  J  r2  sin"  9  dtp2 


=  0 


(5.54) 


The  symmetry  of  the  problem  also  requires  that  the  potential  V  not  vary  with  respect  to  9 
or  0;  in  other  words: 


dv_dv_ 
d9  d(j) 

Thus,  V  is  a  function  of  r  only,  and  Laplace’s  Equation  simplifies  to: 


(5.55) 


hrr(rfrV  M)=° 


(5.56) 


Multiplying  through  by  r2  and  then  integrating  with  respect  to  r,  we  obtain: 

r2^-V(r)  =  C  ,  (5.57) 

dr 

where  C  is  an  arbitrary  constant.  Now  dividing  through  by  r2  and  integrating  with  respect 
to  r  again,  we  obtain: 

V(r)  =  -—  +  C2,  (5.58) 

r 

where  C\  and  C2  are  constants  that  can  be  determined  by  boundary  conditions.  Applying 
the  boundary  conditions,  we  obtain: 

C 

V(r  —  1  m)  = - —  +  C’2  =  100  V  ,  and  (5.59) 

1  m 

V(r  =  3  m)  =  +  C2  =  20  V  .  (5.60) 

3  m 

A  simple  way  to  solve  for  C\  and  C2  here  is  simply  to  subtract  the  second  equation  from 
the  first  equation,  which  eliminates  C2,  then  solve  for  C\  and  use  that  result  to  solve  for  C2. 
One  finds  C\  =  —120  V-m  and  C2  =  —20  V.  Thus: 


V{r) 


120  V  •  m 

r 


-20  V 


1  m  <  r  <  3  m. 


(5.61) 


Note  that  an  answer  like  aV{r)  =  120/r  —  20”  is  dangerously  ambiguous,  unless  you  specify 
as  part  of  the  answer  that  r  must  be  in  meters  and  V  will  be  in  volts. 
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[m0067]  [3] 

5.15-3 


The  symmetry  of  this  problem  suggests  a  solution  in  spherical  coordinates.  Laplace’s  Equa¬ 
tion  in  spherical  coordinates  is: 


1  d 


V2K  =  —  —  r2— -  + 


,dV\ 


r2  dr  \  dr  J  r2  sin  9  dO 


df  dV\ 
sm  9  — —  + 


d2V 


d9  J  r2  sin"  9  dtp2 


=  0 


(5.62) 


The  symmetry  of  the  problem  also  requires  that  the  potential  V  not  vary  with  respect  to  9 
or  0;  in  other  words: 

dV  dV  , 

(5,63) 

Thus,  V  is  a  function  of  r  only,  and  Laplace’s  Equation  simplifies  to: 


hwr  (r2S:*/(r))  =  ° 

Multiplying  through  by  r2  and  then  integrating  with  respect  to  r,  we  obtain: 

r2^)  =  C  , 


(5.64) 


(5.65) 


where  C  is  an  arbitrary  constant.  Now  dividing  through  by  r2  and  integrating  with  respect 
to  r  again,  we  obtain: 

V(r)  =  -—  +  C2,  (5.66) 

r 

where  C\  and  C2  are  constants  that  can  be  determined  by  boundary  conditions.  One  bound¬ 
ary  condition  is  obtained  from  the  surface  of  the  sphere: 


V(r  —  2  m)  = 


Ci 


2  m 


+  Co  =  20  V 


(5.67) 


The  other  boundary  condition  is  obtained  by  noting  that  V(r)  must  go  to  zero  as  r  — »  00. 
Thus: 

V(r  — >  00 )  =  0  +  C2  =  0  (5.68) 

so  C2  =  0  and  C\  =  —40  V-m. 


V(r)  =  4 


40  V  ■  m 


r  >  2  m 


(5.69) 


Note  that  an  answer  like  uV(r)  =  40  jr  V”  is  dangerously  ambiguous,  since  the  units  of  the 
constant  “40”  are  not  clear.  (It  is  OK  -  albeit  tedious  -  to  say  uV(r)  =40 /r  V”  if  you  also 
specify  that  r  is  in  meters.) 
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[m0068]  [1] 

5.16-1 


The  symmetry  of  this  problem  suggests  a  solution  in  cylindrical  coordinates.  Laplace’s 
Equation  in  cylindrical  coordinates  is: 


ld_  (  dV\  &V_ 

p  dp  v  dp  )  p2  d(p2  dz 2 


(5.70) 


The  symmetry  of  the  problem  also  requires  that  the  potential  V  not  vary  with  respect  to  <f> 
or  z ;  in  other  words: 


dv_dv_ 
dcj)  dz 

Thus,  V  is  a  function  of  p  only,  and  Laplace’s  Equation  simplifies  to: 


(5.71) 


d_  (  avA 
dp  V  dp  ) 


0 


(5.72) 


Integrating  both  sides  with  respect  to  p,  we  obtain: 

%v(p)  =  Cl  ■  (5-73) 

where  C  is  an  arbitrary  constant.  Now  dividing  through  by  p  and  integrating  with  respect 
to  p  again,  we  obtain: 

V{p)  =C1\np  +  C2  ,  (5.74) 

where  C\  and  C2  are  constants  that  can  be  determined  by  boundary  conditions.  Applying 
the  boundary  conditions,  we  obtain: 


V(p  —  0.001  m)  =  Ci  In  (0.001  m)  +  C2  =  50  mV  ,  and  (5.75) 

V(p  =  0.002  m)  =  Ci  In  (0.002  m)  +  C2  =  20  mV  .  (5.76) 

A  simple  way  to  solve  for  C\  and  C2  here  is  simply  to  subtract  the  second  equation  from 
the  first  equation,  which  eliminates  C2l  then  solve  for  C\  and  use  that  result  to  solve  for  C2. 
One  finds  C\  =  —43.3  mV  and  C2  =  —249.0  mV.  Thus: 


V(p)  =  -  (43.3  mV)  In  (—)  -  249.0  mV 

VI  m  / 


1  mm  <  p  <  2  mm. 


(5.77) 


Note  that  an  answer  that  does  not  include  “1  m”  in  the  denominator  of  the  argument  of  the 
“In”  function  is  dangerously  ambiguous  unless  you  also  specify  that  p  must  be  in  meters. 
Here  is  another  correct  solution,  this  time  with  the  argument  in  units  of  millimeters: 


V(p)  =  —  (43.3  mV)  In  ( — — — +  50.0  mV  ,  1  mm  <  p  <  2  mm.  (5.78) 

VI  mm/ 

Note  that  C2  depends  on  the  units  of  p  in  the  argument  of  the  logarithm  function. 
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[m0021]  [1] 

5.18-1 


A  good  way  to  get  the  charge  density  is  to  first  find  the  electric  held,  and  then  to  apply 
the  boundary  condition  that  relates  electric  held  to  surface  charge  density  on  a  conducting 
surface.  The  electric  held  intensity  is,  in  general: 


E  =  -W 


~ldV  _ „dV 
p  d(j)  Z  dz 


(5.79) 


The  last  two  terms  are  zero  because  the  the  answer  cannot  vary  with  respect  to  <ft  or  z.  So 
we  have: 


E 


(43.3  mV)  In 


-  249.0  mV 


„  43.3  mV 
P  - 


P 


(5.80) 


The  relevant  boundary  condition  on  the  inner  conductor  is  that  the  normal  component  of  the 
electric  hux  density  D  equals  the  surface  charge  density.  The  normal  to  the  inner  conductor 
is  +p,  so  we  have: 


Ps  —  (+p)  '  D|p=i  mm 


P  ■  eE|p=i  mm  =  ere0  (43.3  V/m)  . 


(5.81) 


Since  e0  =  8.854  x 
inner  conductor  is 


10  12  F /m  and  er 
+804  pC/m2  . 


2.1,  we  have  that  the  surface  charge  density  on  the 
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[m0021]  [2] 

5.18-2 


(a)  Summarizing  the  problem  statement,  we  have  a  sphere  of  radius  a  —  2  m  containing 
uniformly- distributed  charge  with  volume  density  pv  =  3  pC/m3,  and  the  media  is  a  dielec¬ 
tric  with  er  =  4.5  everywhere.  Poisson’s  Equation  is 

VV  =  ~  (5.82) 

Note  that  you  could  also  do  this  problem  by  integrating  over  the  charge  distribution,  and 
that’s  a  great  check.  However,  the  problem  statement  requires  you  to  use  Poisson’s  Equation. 
The  symmetry  of  the  problem  suggests  the  use  of  spherical  coordinates.  Noting  that  dV/d9 
and  dV/dej)  should  be  zero  due  to  symmetry,  we  find 


It’s  straightforward  to  solve  for  V 


where  C\  is  an  arbitrary  constant. 


1  d 

7^~dr  '  r 


,dV_ 

dr 


Pv 

e 


in  this  case.  Here  we  go: 
d 


dr 


2dV_\  _ 
T  dr  ) 


Pv  2 
— r 
e 


,dV 


2  “  '  Pv  3  . 

r  —  =  r*  +  C  i 
dr  3e 

Continuing: 


9V  =  _Pvr,C i 

dr  3e  r2 


(5.83) 


(5.84) 

(5.85) 


(5.86) 


V(r)  =  -^r2  +  ^  +  C3  (5.87) 

6e  r 

where  C2  and  C3  are  arbitrary  constants.  At  this  point  you  should  confirm  this  result  by 
making  sure  it’s  a  solution  to  the  original  equation,  and  also  by  checking  units. 

Outside  the  sphere,  pv  =  0.  Thus: 


D 

V(r)  —  —  +  A0  r>a  (5.88) 

r 

Here  we  have  replaced  the  constants  C-2  and  C3  with  Ba  and  Aa  respectively.  This  is  to  remind 
us  that  the  constants  may  be  different  should  we  consider  the  region  inside  the  sphere  (as 
we  shall  soon  do).  We  can  determine  the  value  of  the  constant  A0  by  noting  that  V(r)  — >  0 
as  r  — »  oo,  since  the  total  charge  is  finite  and  contained  within  a  finite  region.  Therefore, 
Aa  must  be  zero,  leading  to 

V(r)  —  —  r>a  (5.89) 

r 

To  determine  the  value  of  the  constant  Ba  we’re  going  to  have  to  make  some  kind  of  con¬ 
nection  with  V (r)  inside  the  sphere.  Inside  the  sphere: 


V(r)  = 


77-r2  +  —  +  Ai  r  <  a 
be  r 


(5.90) 
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Here  we  have  replaced  the  constants  C2  and  C3  with  Bi  and  Ai  respectively.  We  can  deter¬ 
mine  the  value  of  the  constant  B {  by  noting  that  V(r)  must  be  finite  as  r  — y  0,  since  the 
charge  density  is  finite  at  r  =  0.  Therefore,  B j  must  be  zero,  leading  to 

V(r)  =  — — r2  +  Ai  r<a  (5.91) 

6e 

Now  we  apply  the  boundary  condition  at  the  surface  of  the  sphere.  Note  that  there  is  no 
requirement  for  potential  to  be  continuous  (and  it  wouldn’t  do  us  any  good  even  if  there 
were,  since  we’d  be  stuck  with  one  equation  and  two  unknowns).  The  relevant  boundary 
condition  at  r  =  a  is  that  the  normal  component  of  the  electric  field  should  be  continuous: 
Specifically, 

[D0(r  =  a)  -  Dj(r  =  a)]  •  r  =  ps  (5.92) 

where  D0  and  D,  are  the  electric  flux  densities  outside  and  inside  the  sphere  respectively, 
and  ps  is  the  surface  charge  density.  The  surface  charge  density  ps  is  zero,  since  all  the  charge 
is  taken  into  account  as  the  volume  charge  density  pv.  Also,  D  =  eE;  therefore,  we  have 

[E 0{r  =  a)  —  Ej(r  =  a)]  ■  r  =  0  (5.93) 


Next  we  note  E  everywhere  should  be  oriented  in  the  r  direction  due  to  symmetry.  Thus, 
we  find: 

E0{r  =  a)  =  Ei(r  =  a)  (5.94) 

We  can  find  the  electric  flux  density  by  taking  the  gradient  of  the  potential: 


E  =  -W  =  -r^ 
or 


—S  +  A 

6e 


~  Pv  , 

=  r — r  r  <  a 
3e 


E  =  -W  =  -r^- 
or 


Bn 


Bn 


=  r- 


r  >  a 


Now  applying  the  boundary  condition  ( Equation  15.941): 


(5.95) 

(5.96) 


!“=§  <5J7> 

Solving  for  B0  and  substituting  the  result  back  into  Equation  15.891  we  obtain: 

V(r)  =  - - =  — -  r>a  (5.98) 

3er  3ereo r 

(Good  time  for  a  units  check!)  Finally,  the  answer: 

,  .  0.201  V-m 

V(r)  = -  r>a  (5.99) 


(b)  From  the  previous  equation,  H(3  m)  =  66.9  mV. 
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[m0021]  [3] 

5.18-3 


A  good  way  to  get  the  charge  density  is  to  first  find  the  electric  held,  and  then  to  apply 
the  boundary  condition  that  relates  electric  held  to  surface  charge  density  on  a  conducting 
surface.  The  electric  held  intensity  is: 


X.  ___  t  r  *dV  AdV  ,  1  dV 

E  =  _W  =  -rA - 6~  AT  _  ^ — —qAa  ■ 

or  r  at)  r  sm  9  ocp 


(5.100) 


The  last  two  terms  are  zero  because  the  the  answer  cannot  vary  with  respect  to  9  or  </>.  So 
we  have: 

„  A dV  „  d 
~  ~ T~d r  ~ 


120  V  •  m 


20  V 


=  +r: 


120  V  •  m 


(5.101) 


The  relevant  boundary  condition  on  the  inner  conductor  is  that  the  normal  component  of  the 
electric  hux  density  D  equals  the  surface  charge  density.  The  normal  to  the  inner  conductor 
is  +r,  so  we  have: 


ps  =  (+?)  •  D|r=i  m  =  f  ■  e0E|r=i  m  =  e0  (120  V/m)  .  (5.102) 

Since  eo  =  8.854  x  10^12  F/m  in  free  space,  we  have  that  the  surface  charge  density  on  the 
inner  conductor  is 


1.06  nC/m2 
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[m0112]  [1] 

5.22-1 


The  net  charge  in  the  capacitor  is  zero  because  the  charges  on  the  two  plates  is  equal 
and  opposite.  The  charge  on  the  positively-charged  plate  is 


Q+ 


CV  =  (20  pF)  (3  V) 


+60  pC 


84 


[m0070]  [1] 

5.23-1 


a)  The  equivalent  circuit  is  just  a  resistor  R  in  parallel  with  a  capacitor  to  which  we’ll 


assign  the  variable  Cx.  The  impedance  of  the  capacitor  is  — j/uCx .  Thus,  the  magnitude 
of  the  capacitor’s  impedance  decreases  with  increasing  frequency.  The  total  impedance  is 
the  parallel  combination  of  R  =  200  and  Cx.  The  effective  resistance  will  decrease  with 
increasing  frequency. 


(b)  This  structure  looks  a  lot  like  a  parallel  plate  capacitor.  Neglecting  fringing  fields, 
capacitance  is  estimated  as 


a  = 


eHW  6q erHW  (8.854  x  10  12  F/m)  •  37  •  (0.3  mm)  (0.3  mm) 


0.6  mm 


49.1  fF 


(c)  The  impedance  is  R  ||  (—j/uCx).  At  /  =  10  GHz,  we  have  144.8  —  j89.4  fl,  so  the 


effective  resistance  is  144.8  Q  . 
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[m0070]  [2] 

5.23-2 


From  the  problem  statement:  C  <  3  pF,  d  =  2  mm,  and  er  =  3.0.  Note: 

eA  £r£oA 

where  A  is  the  area  in  common.  So: 

A  <  (3  PF)^ 
ereo 

Therefore,  the  common  area  must  be  <  2.26  x  10~4  nr. 
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[m0070]  [3] 

5.23-3 


(a)  Given  e(z) 
unit  less  . 


e0  (az  +  b)  is  permittivity,  the  units  of  a  must  be 


and  b  must  be 


(b)  This  problem  is  really  quite  similar  to  the  derivation  presented  in  the  book.  As  in  that 
derivation,  electric  flux  density  between  the  plates  is 

D  «  -z pS)+  (5.103) 

where  pS)+  is  the  charge  density  on  the  positively-charged  plate  at  z  —  d.  The  electric  flux 
density  is  unchanged  in  this  problem  because  electric  flux  does  not  depend  on  the  material 
in  which  it  exists  (unlike  the  electric  field  intensity,  E).  When  it  comes  time  to  compute  the 
potential  across  the  plates,  we  find: 


G  = 


E  •  dl 


I  z= 0 

f 

lz= 0 


D 


-n  \  eiz 


Ps 


■  (z dz) 

~  Ps,+ 

_ 2 _ - _ 

£q  | az  +  b | 

dz 


■  (z  dz) 


e0  Jz= o  az  +  b 

You  can  solve  the  integral,  or  just  look  it  up  in  a  table  of  integrals.  Continuing: 


V  = 


Ps,+ 


eo  \  a, 

Ps,+  1 
eo  a 
Ps,+  1 
e0  a 
Ps,+  1 
e0  a 


(5.104) 

(5.105) 

(5.106) 

(5.107) 


-  In  \az  +  b\ 

(5.108) 

a 

0/ 

(In  \ad  +  b\  - 

In  \b\) 

(5.109) 

(5.110) 

Kf+i) 

(5.111) 

We  can  dispense  with  the  absolute  value  operator  above  since  the  argument  is  always  non¬ 
negative.  Finally: 

Ps,+A 


A  Q-\ 


C  =  — 


V  (ps,+/eo)  (1/a)  111  ( ad/b  +  1) 


which  simplifies  to: 


C  &  e0A 


In  ( ad/b  +  1) 


(5.112) 

(5.113) 
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(c)  Units  check:  eo  has  units  of  F/m,  A  has  units  of  m2,  and  a  has  units  of  1/m.  The  “In 
factor  in  the  denominator  is  unitless.  Thus  we  find  that  C  has  units  of  F,  as  expected. 


(d)  Recall  that  e(z)  =  e0(az  +  b).  If  the  permittivity  is  uniform,  then  a  must  be  zero. 
Subsequently  b  must  be  the  relative  permittivity,  er.  Equation  15. 1 131  becomes: 


C  ~  e0A  lim 


n->o  ln(ad/er  +  l) 


(5.114) 


Note  we  have  to  be  careful  because  both  numerator  and  denominator  are  going  to  zero. 
Applying  L’Hopital’s  Rule,  we  take  the  derivative  with  respect  to  a  of  the  numerator  and 
denominator: 


lim 


=  lim 


o  ln(ad/er  +  l)  a-+o  (d/er)  /  (ad/er  +  1)  d 


(5.115) 


Substituting  this  result  into  Equation  15.1141  we  obtain: 


C^e0A 


6r 

~d 


eA 

~d 


(5.116) 


which  is  the  expected  result  (i.e.,  the  one  we  had  already  derived  for  uniform  permittivity). 


(e)  In  terms  of  the  variables  established,  A  =  400  pm2,  d  =  0.5  mm,  5  =  2,  and  a  = 
(10  —  2 )/d  =  16000  m_1.  (You  should  check  that  this  choice  of  a  and  b  gives  you  e  =  2eo  at 
z  —  0  and  e  =  10eo  at  z  —  d.)  Equation  15. 1 131  becomes: 


C  «  e0A 


a 

In  (5) 


35.2  pF 


(5.117) 


[m0113]  [1] 

5.24-1 


From  the  problem  statement,  C'  =  30  pF/m  and  er  =  2.25  for  polyethylene.  The  capacitance 
of  the  original  coaxial  cable  is 

,  _  27res 
In  (b/ a) 

where  es  =  ere0  for  the  spacer  material  and  b/a  is  the  ratio  of  the  radius  of  the  outer 
conductor  to  that  of  the  inner  conductor.  The  capacitance  of  a  coaxial  cable  that  is  identical 
except  polyethylene  is  replaced  with  air  (es  =  eo)  is: 


C 

w  new 


27T6Q 

In  (b/a) 


Comparing  the  two  equations,  we  observe: 


C 

w  new 


C' 

6r 


13.3  pF/m 
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[m0113]  [2] 

5.24-2 


From  the  problem  statement,  a  =  1  mm,  6  =  3  mm,  es  =  eo,  and  V  =  +1.5  kV  mea¬ 
sured  at  the  outer  conductor  relative  to  the  inner  conductor.  The  capacitance  of  this  cable 

is 


C' 


27Te0 
In  (6/a) 


=  50.6  pF /m 


The  outer  conductor  is  positively-charged,  and  the  line  charge  density  on  this  conductor  is 


Pi 


C'V  + 


+76.0  nC/m 


The  circumference  of  the  outer  conductor  is  27t6,  so  the  surface  charge  density  is 


Ps 


Pi 

2nb 


+4.03  /iC/m2 
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[m0114]  [1] 

5.25-1 


From  the  problem  statement,  C  =  4.7  mF  and  V  =  16  V.  The  energy  stored  is 

1 


Wp  =  -CV1  ^  602  mJ 
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[m0114]  [2] 

5.25-2 


From  the  problem  statement,  C 
ideal  parallel  plate  capacitor, 


3.5  pF,  d  =  0.1  mm,  er 


C  = 


eA 

~d 


where  e  =  ereo  is  the  spacer  permittivity  and  A  is  the  plate 
we  find: 

A  = - =  3.95  x  10'6  m2 


£o£r 


The  volume  of  the  capacitor  is 


Ad  =  3.95  x  10'10  m3 


The  energy  in  the  capacitor  is 

We  =  X-CV2  =  15.75  pJ 

Therefore,  the  energy  density  is 


We 

Ad 


39.8  mJ/m3 


=  10,  and  V  =  3  V.  For  an 


area.  In  the  present  problem, 
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Chapter  6 

Steady  Current  and  Conductivity 
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[m0071]  [1] 

6.4-1 


(a)  Resistance  Rsteei  Per  unit  length  l  of  the  steel-only  wire: 

t-,/  Rsteei  1  1 

-t^steel 


l  @ steel  •  ^a?  (1.00  x  106  S/m)  ■  n  (0.1  mm)' 


31.8  O/m 


(b)  Resistance  per  unit  length  of  gold  clad  having  outer  radius  b : 

/?'  1 

9°ld  agold-n  (ft2 -a2) 


The  total  resistance  per  unit  length  R'totai  =  10  O/m  is  the  parallel  combination: 


1 


total 


l  i  l 

■n' steel  gold 


1 


RL 


steel 


+  Vgoid  -7 :  (b2  -  a2) 


Solving  for  b: 


b  = 


p/ 

_  ^ total 


R' 


T  C 7 gold  ' 


steel 


® gold 


=  0.10263  mm 


•  7 r 


So  the  required  thickness  of  gold  is  b  —  a,  which  is  2.63  //m 
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[m0071]  [2] 

6.4-2 


In  the  transmission  line  equivalent  circuit  {R! ,G' ,C ,L')  model,  R'  is  a  series  resistance.  Also, 
any  current  applied  to  either  conductor  must  return  on  the  other  conductor.  Therefore,  we 
have 

R1  =  Kc  +  R'oc  (6-i) 

where  R'ic  is  the  resistance  per  length  of  the  inner  conductor  and  R'oc  is  the  resistance  per 
length  of  the  outer  conductor.  Note 


where  a ic  is  the  inner  conductor  conductivity  and  Aic  is  the  cross-sectional  area  of  the  inner 
conductor.  Thus,  R'ic  =  0.164  O/rri.  Also 


R'  = 


&nr.An 


(6.3) 


where  aoc  =  (JjC  (from  the  problem  statement)  and  Aoc  is  the  cross-sectional  area  of  the  outer 
conductor,  through  which  the  current  flows.  Note: 


Aoc  =  ■Kb 2  -  Kb\ 


(6.4) 


where  b\  and  62  are  the  radii  of  the  inner  and  outer  surfaces,  respectively,  of  the  outer 
conductor.  From  the  problem  statement  we  have 

b\  =  ^1  —  ^  —  0.1809  cm  (6.5) 

b2=  (l  +  b  =  0.1901  cm  (6.6) 

(Check:  the  mean  of  bi  and  b2  is  (fex  +  b2)  /2  =  b,  as  expected.)  So  Aoc  =  1.081  x  10~6  m2, 
and  subsequently,  R'oc  =0.0406  O/m.  Finally,  we  obtain  R'  =  R'ic  +  R'oc  =  0.205  O/m. 
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[m0071]  [3] 

6.4-3 


If  the  voltage  drop  is  to  be  reduced  by  a  factor  of  2,  then  the  resistance  must  be  decreased 
by  a  factor  of  2.  The  DC  resistance  of  a  wire  is  l /a A  where  l  is  length,  a  is  conductivity, 
and  A  is  cross-sectional  area.  The  use  of  the  term  “diameter”  implies  the  wire  has  circular 
cross  section,  so  the  original  wire  has  A  =  n  (D0/ 2) 2  and 


<77r(T>o/2)2 


For  this  to  be  reduced  by  a  factor  of  two,  the  new 


diameter  must  be 


V2-D0 
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[m0071]  [4] 

6.4-4 


The  DC  resistance  of  such  a  resistor  is  R  =  l /a  A  where  A  is  cross-sectional  area.  In  this 
case,  we  may  write 

R  1  l  ,4/ 

crA  cm  (D/2)2  cm 

where  D  is  diameter.  Note  that  D  should  increase  by  a  factor  of  y/2  in  order  to  reduce  R 
by  a  factor  of  two.  Thus,  D  becomes  =  1.41  mm. 
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[m0105]  [1] 

6.5-1 


The  conductance  per  unit  length  is 


G' 


‘I'KOg 

In  {b/ a) 


where  as  is  the  spacer  conductivity,  and  a  and  b  are  the  radii  of  the  inner  and  outer 
conductors,  respectively.  From  the  book,  RG-59  has  <js  =  5.9  x  10-5  S/m  and  exhibits 
G'  =  200  /j, S/m  normally.  From  the  appendix  “Conductivity  of  Some  Common  Materials,” 
crs  ~  5  S/m.  The  worst  case  is  that  spacer  assumes  the  much  higher  conductivity  of  seawater, 
in  which  case: 

G'  (200  /i S/m)  5  9  / g/m  = 


17.0  S/m 
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[m0105]  [2] 

6.5-2 


Let  the  ground  plane  be  at  z  =  0,  and  let  the  trace  be  at  z  =  h.  From  the  problem 
statement,  the  current  /  is  positive  when  flowing  into  the  trace  and  from  the  ground  plane. 

(a)  Under  the  condition  that  W  h  we  assume  that  most  of  the  current  in  the  transmission 
line  flows  directly  from  the  trace  to  ground  plane  in  the  — z  direction,  and  that  the  fraction 
of  current  that  does  not  satisfy  this  condition  (i.e. ,  the  current  close  to  the  edges  of  the 
trace)  is  negligible.  Thus,  we  are  justified  in  assuming  the  current  density  is  approximately 
uniform  throughout  the  region  directly  underneath  the  trace.  Therefore  the  magnitude  of 
the  current  density  is  approximately  total  current  /  divided  by  trace  area  Wl ,  where  l  is  the 
length  of  the  trace.  Under  this  same  approximation,  the  magnitude  of  the  current  density 
is  assumed  to  be  zero  beyond  the  trace.  Summarizing: 


f-z  I/Wl, 

directly  underneath  trace;  and 

J~i  o, 

otherwise. 

(6.7) 


where  l  is  the  length  of  the  trace. 


(b)  The  electric  field  intensity  is  given  by  Ohm’s  law: 

J  I  — z I/Wlas,  directly  underneath  trace;  and 

crs  I  0,  otherwise. 


Subsequently: 


V  = 


(6.8) 


(6.9) 

(6.10) 


Note  the  start  point  is  z  =  0,  since  this  is  the  negative  terminal  with  respect  to  a  current 
source  driving  the  transmission  line.  Similarly  the  end  point  is  z  =  h,  since  this  is  the 
positive  terminal  with  respect  to  a  current  source  driving  the  transmission  line.  Continuing: 


V  « 


I 

Whj~s 


(6.11) 


The  integral  is  equal  to  h.  Thus: 


V  « 


Ih 

Wlas 


(6.12) 


This  is  the  potential  measured  at  the  trace  relative  to  the  potential  at  the  ground  plane. 
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(c)  Conductance  G  is  determined  as  follows: 


G  = 


Wl<7s 

h 


Conductance  per  unit  length  is  G'  =  G /l,  so 


(6.13) 


(6.14) 


(d)  Since  the  trace  and  ground  plane  are  specified  to  be  perfectly  conducting,  R!  =  0  and  the 
only  physical  mechanisms  to  consider  are  G',  L',  and  C' .  At  DC,  There  is  no  contribution 
from  L'  since  it  is  in  series  with  the  trace,  and  there  is  no  contribution  from  C  since  it 
connects  trace  to  ground  plane.  Therefore  Z  =  V/ 1  =  1/G.  From  Equation  16.131  we  obtain: 


(6.15) 
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[m0106]  [1] 

6.6-1 


From  the  problem  statement: 

Length  l  =  1.2  cm, 
radius  a  =  1.6  mm, 

J  is  uniform  (constant)  in  the  resistor, 

E  =  zEq / \J~p  where  E0  =  3  V  ■  m-1/2,  and 
P  =  5  W. 

Let  a  be  the  conductivity  of  the  material  comprising  the  resistor.  Then 

P  =  [  a  |E|2  dv 
Jv 

where  V  is  the  volume  representing  the  resistor. 


We  cannot  assume  the  material  comprising  the  resistor  is  homogeneous.  So,  what  do  we 
know  about  a?  Recall  Ohm’s  Law,  J  =  crE.  Since  E  is  proportional  to  1  j  ^fp  and  J  is  inde¬ 
pendent  of  p,  a  must  have  the  form  cr  =  (X0v/p  where  <70  is  a  constant  having  units  of  S-m-3/2. 


Continuing, 

P  =  I  (<7oy/p)  dv  =  a°E0  J  P~l/2dv 

Let  us  assume  the  ends  of  the  resistor  are  at  z  =  0  and  z  —  l.  Then: 

pa  p2n  pi 

P  =  a0E2  /  /  /  P~lf 2  [dp  ( pd(p )  dz] 

J  p= 0  J  <j)= 0  J  z= 0 


Solving  for  cr0: 


and  subsequently, 


^o-^o 


n2n 


/  P+1/2dp)  [  I  dcj) 

Ip= 0  J  \Jcf>=0 


=  OqEI  0a3/2^)  (2?r)  (/) 


dz 


<z= o 


3  P 


(T0  = 


4vra3/2/E2 


^  173  kS  ■  m~3/2 


(173  kS  ■  m  3/2)v^ 
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Chapter  7 
Magnetostatics 
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[m0115]  [1] 

7.1-1 


Divergence  of  the  electric  field: 

V  •  D  =  pv  ,  so 

V  ■  eE  =  pv 

where  e  has  units  of  F/m  and  pv  has  units  of  C/m3. 
Curl  of  the  electric  field: 

V  x  E  =  0 


Divergence  of  the  magnetic  field: 

V  •  B  =  0  ,  so 

V  ■  H  =  0 

Curl  of  the  magnetic  field: 

V  x  H  =  J 


Noting  J  =  crE: 

V  x  H  =  ffE 

where  a  has  units  of  S/m. 


103 


[m0047]  [1] 

7.3-1 


The  differential  form  of  Gauss’  Law  for  magnetism  requires 

V-  B  =  0 

According  the  reported  measurement: 

(r\  r\  r\  \ 

Xt-  +  +  Ztt  )  •  ZB0x2  =  2 B0x 

ox  ay  oz ) 


Therefore,  the  measurement  is  plausible  only  if  B0  is  zero. 
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[m0119]  [1] 

7.5-1 


The  magnetic  flux  density  (not  the  same  as  magnetic  flux!)  in  this  case  is  given  by 

B(P)=^.  (7.1) 

The  magnetic  flux  is  simply  B  integrated  over  the  area  of  the  loop  (i.e.,  flux  divided  by  area, 
times  area,  is  flux): 


/»  /*23  cm  /*30  cm  T 

$  =  /  B  •  ds  =  /  /  ■  4>  dp  dz  (7.2) 

J S  J  p=3  cm  J  Z=0  ‘i'K p 

Where  S  is  the  area  enclosed  by  the  loop,  and  the  absolute  values  of  z  don’t  matter  due  to 
symmetry,  so  you  can  pick  any  z' s  you  want  as  long  as  you  cover  30  cm  in  the  z  direction. 
Evaluating: 


$  = 


IM)I 

2tt 


IJ  p=3 


23  cm  dp 

P  . 


Solving  for  /  we  have 

/  =  (3  pT-m2) 


-30 


Uz= 0 


2tt 


dz 


47t  x  10-7  H/m 


Pd 

2v r 


In® 

3 


in  — 
3 


(30  cm)  =  3  /xT-m2 


(7.3) 


-l 


.  _]  T-m" 

(0.3  m)  1  =  24.5  — — 

H 


(7.4) 


To  get  to  units  of  A,  the  traditional  units  of  current,  recall  that  inductance  (H)  is  defined 
as  magnetic  flux  (T-m2)  divided  by  current  (A),  so  we’re  already  there!  Thus,  we  have 


1=  24.5  A 
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[m0119]  [2] 

7.5-2 


The  magnetic  flux  density  (not  the  same  as  magnetic  flux!)  is  given  by 


=  <7-5> 

The  magnetic  flux  is  simply  B  integrated  over  the  area  of  the  loop  (i.e.,  flux  divided  by  area, 
times  area,  is  flux): 


r-0.02 


r-0.10 


<f>  =  /  B-ds  = 


I  P=0.01  m  J  Z= 0 


■  IM)I 

2np 


(j)  dp  dz 


Evaluating: 


$  = 


p0I 


r-0.02 


dp 


2tt  Jp= o.oi  m  P 


r-0.10 


lJz= 0 


dz 


(4vr  x  10"7  H / m)  (3  A) 


2tt 


In 


0.02 

OOl 


(0.10  m)  =  41.6  nT-m2 


Note  that  this  may  also  be  written  as  41.6  nWb  ,  since  1  Wb  =  1  T/m2. 


(7.6) 


(7.7) 

(7.8) 


106 


[m0119]  [3] 

7.5-3 


Ampere’s  Law  is 


H  •  dl  =  I, 


encl 


(7.9) 


where  C  is  any  path  which  encloses  the  current.  A  convenient  path  is  just  a  constant-;?  circle 


with  radius  a: 


r27 r 

/  H  •  (j)  a  d(j)  =  I 

Jtj)=0 

(7.10) 

Also,  H  =  B / fiQ,  so 

/*2n 

/  (J)Jqci  ■  4> 

J<j>= 0 

1 — i 

o 

S3 

(7.11) 

Thus: 

/  =  27ra2J0  = 

20.0  mA 

(7.12) 
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[m0119]  [4] 

7.5-4 


Apply  the  right  hand  rule.  When  the  thumb  of  the  right  hand  points  in  the  —  y  direc¬ 
tion,  the  curled  fingers  of  the  right  hand  point  in  the  +z  direction  at  (+1,  +1,0)  m. 
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[m0119]  [5] 

7.5-5 


The  measurements  are  explained  by  a  wire  aligned  along  the  y- axis,  with  current  flowing 
in  the  —  y  direction.  This  can  be  confirmed  using  the  right-hand  rule  -  point  the  thumb  of 
your  right  hand  in  the  direction  of  the  current,  and  the  curled  fingers  of  your  right  hand 
point  in  the  direction  of  the  magnetic  field. 
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[m0120]  [1] 

7.6-1 


From  the  problem  statement,  we  have  that  Coil  1  has  Ni  =  100  and  I\  —  2  A.  Coil  2 
has  N\  =  300  and  I2  =  4  A,  and  is  wound  in  the  opposite  direction.  Both  coils  have 
l  =  10  cm  and  ji  =  /io-  For  Coil  1  we  have 


Bi  —  bi/io 


Nj! 

I 


where  bx  points  in  the  direction  of  Bx  inside  the  coil.  For  Coil  2  we  have 

N2I2 


B2  —  — bi/iQ- 


The  total  field  B  =  Bi  +  B2,  so: 


t>i  •  (Bi  +  B2 


/^O" 


N2I2 


A4  o- 


l 


A4  0 


=  y  I Nih  -  N2I2 1  =  12.6  mT 
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[m0049]  [1] 

7.7-1 


The  magnetic  field  in  either  a  straight  coil  or  a  toroidal  coil  is  proportional  to  current. 
Current  is  proportional  to  the  conductivity  of  the  wire  forming  the  coil.  Therefore,  doubling 
the  conductivity  will  double  the  magnetic  field  strength. 
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[m0121]  [1] 

7.8-1 


(a)  Positive  V  corresponds  to  current  flowing  in  the  +z  direction  along  the  trace.  The 
direction  of  the  associated  magnetic  field  can  be  determined  using  the  following  “right  hand 
rule:”  Orient  the  thumb  of  your  right  hand  in  the  reference  direction  of  current  flow  in  the 
trace,  and  observe  the  direction  in  which  the  curled  fingers  of  your  right  hand  point.  We  see 
that  deep  inside  the  transmission  line,  the  direction  is  +x  . 


(b)  The  integral  form  of  Ampere’s  law  is: 

H  •  dl  =  Iencl  (7.13) 

The  only  “hard”  requirement  on  C  is  that  it  enclose  some  of  the  relevant  current.  The  most 
convenient  choice  for  C  is  shown  below: 


X 


Here  are  the  considerations  leading  to  this  choice: 

•  The  indicated  direction  of  C  is  consistent  with  the  expected  direction  of  the  magnetic 
field,  as  determined  in  part  (a). 

•  We  choose  a  path  that  lies  in  plane  of  constant  z,  since  this  minimizes  the  number 
of  varying  parameters  required  to  describe  the  path.  The  precise  choice  of  z  is  not 
important  as  long  as  it  is  as  far  from  either  end  of  the  transmission  line,  where  we 
would  expect  fringing  fields  to  become  potentially  important. 

•  Segment  A  is  a  line  of  constant  y  which  lies  entirely  within  the  transmission  line  (i.e. , 
between  y  —  0  and  y  =  h),  and  is  required  since  we  need  some  portion  of  C  to  be 
coincident  with  the  location  where  we  wish  to  determine  the  field. 
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•  Segment  C  is  chosen  to  lie  along  a  line  of  constant  y  at  oo.  This  is  convenient  because  we 
expect  the  magnetic  field  go  to  zero  as  the  distance  from  this  finite  structure  increases 
to  infinity. 


•  Segments  B  and  D  are  chosen  to  lie  along  lines  of  constant  x  so  as  to  connect  Segments  A 
and  C  using  paths  that  can  be  described  in  the  minimum  number  of  varying  parameters: 
For  these  segments,  the  only  variation  is  along  y.  We  choose  x  =  —W/ 2  and  x  =  +W / 2 
since  this  closes  path  C  with  the  shortest  total  path  length  that  encloses  all  of  the 
current  of  interest.  Choosing  a  path  wider  than  the  trace  would  result  in  integration 
over  a  region  where  more  field  lines  are  significantly  curved. 

(c)  Continuing  with  the  left  side  of  Equation  17. 131  using  the  path  determined  in  part  (b): 


U  dl=  /  H  ■  dl 


H  •  dl 


Hdl  +  /  H  dl 


IB 


1C 


Id 


«  /  H-dl +  0  +  0  +  0  (7.14) 

J  A 

The  integral  over  Segment  C  is  exactly  zero  because  H  =  0  along  this  segment,  as  explained 
in  part  (b).  The  integral  over  Segments  B  and  D  is  approximately  zero  because  H  ■  dl  ~  0 
along  these  segments.  Along  Segment  A,  we  have  H  =  +x.H(x,y).  Also,  Ienci  =  V / R.  Thus 
Equation  17. 131  reduces  to: 


xH(x,  y)}  ■  l+'kdx] 


H{x,y )  •  W 


V 
R 

V 
R 


(7.15) 

(7.16) 


Equation  17.161  indicates  that  the  magnetic  field  along  Segment  A  does  not  depend  on  x  or 
y\  at  least  given  the  assumptions  made  to  this  point.  Therefore  H(x,y)  is  considered  a 
constant.  Thus  we  find  that 


H 


V 

Wr 


deep  inside  the  transmission  line 


(7.17) 


(d)  The  magnetic  flux  density  B  =  pH.  Since  the  spacer  material  is  non-magnetic,  /i 
Thus: 


B 


*WR 


deep  inside  the  transmission  line 


«  +o- 


(7.18) 


(e)  The  differential  form  of  Gauss’  law  for  magnetism  is  V  •  B  =  0.  Divergence  is  essentially 
the  first  derivative  with  respect  to  position.  From  Equation  17. 181  we  see  that  B  is  constant  (at 
least  approximately)  with  position.  Therefore  Equationl7.18lis  consistent  with  the  differential 
form  of  Gauss’  law  for  magnetism.  The  differential  form  of  Ampere’s  law  is  V  x  H  =  J. 
In  the  dielectric  spacer,  J  =  0.  Curl  is  also  essentially  the  first  derivative  with  respect  to 
position.  Thus,  Equation  17. 181  is  consistent  with  the  differential  form  of  Ampere’s  law. 

(f)  Using  Equation  17. 181 


„  (47t  x  10  '  H/m)  (+5  mV) 
(6  mm)  (50  f2) 


x20.9  nT 


(7.19) 
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[m0123]  [1] 

7.12-1 


The  inductance  of  a  linear  inductor  depends  only  on  geometry  and  materials;  therefore, 
the  inductance  remains  1  H  . 
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[m0124]  [1] 

7.13-1 


If  the  loops  are  close  together,  then  presumably  the  magnetic  flux  $  through  each  wind¬ 
ing  is  equal.  Thus, 


A  ^ 

I 


where  N  is  number  of  linkages,  which  is  this  case  is  the  number  of  windings.  Thus, 


$ 


LI  _  LI 
~N  ~  ~2 
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[m0124]  [2] 

7.13-2 


From  the  problem  statement: 
length  l  =  5  cm, 
radius  a  =  5/2  =  2.5  mm, 
number  of  windings  N  =  300,  and 
relative  permeability  /ir  =  200. 


Since  l  a  and  the  winding  density  N/l  is  large,  we  may  us  the  “long  straight  coil” 
expression 

_  /j,N2A 


In  the  present  problem: 


(Hrl^o)  N2  (tt a2) 

l 


8.88  mH 


where  /j,q  =  47T  x  10  7  H/m. 
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[m0127]  [1] 

7.15-1 

Since 


We  have 


We  = 


2  ■  2  mJ 
47  mH 


292  mA 
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[m0127]  [2] 

7.15-2 


The  energy  initially  stored  in  the  inductor  is 

1 


W,m  =  -Lr  =  6  nJ 


The  energy  stored  in  the  capacitor  after  the  transfer  is 

WP  =  lev2  =  6  nJ 


so  V 


1.73  V 
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Chapter  8 

Time- Varying  Fields 
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[m0055]  [1] 

8.3-1 


Faraday’s  Law  says  that  the  potential  (or  “emf”)  induced  in  the  coil  is 

vcm,  =  -Njtm 


where 


$(t)  =  /  B  •  ds 


and  where  B (t)  is  the  magnetic  flux  density  and  S  is  the  surface  defined  by  the  cross-section 
of  the  coil.  When  the  magnetic  flux  is  not  varying  with  time,  the  potential  is  zero  ,  so  this 
is  the  answer  to  parts  (a)  and  (c). 

While  the  magnetic  field  is  being  reduced,  a  non-zero  potential  is  possible.  Since  the  magnetic 
field  is  spatially-uniform  and  parallel  to  the  axis  of  the  coil,  the  above  integral  simplies  to 

$(t)  =  B{t)  ■  A 

where  B(t)  is  the  scalar  magnetic  flux  density  and  A  is  the  cross-sectional  area  of  the  coil. 
We  do  not  know  precisely  how  B  varies  with  time  (i.e. ,  linearly  with  time?  exponential 
decay?  etc.),  so  we  cannot  take  a  formal  derivative.  We  can  however  estimate  the  derivative: 


d  ^ ,  . 


A$  B(to  +  200  ms)  •  A  —  B(to)  ■  A 
At  200  ms 


where  to  is  the  time  at  which  the  magnitude  of  the  magnetic  field  begins  to  decrease.  We 
also  note 

B(t0)  =  nriioH (to)  =  (2  x  105)  •  (47t  x  10”'  H/m)  •  (20.0  niA/m)  =  5.03  mT 
and  B{to  +  200  ms)  =  1.01  mT.  Since  A  =  200  cm2  =  0.020  m2,  we  may  now  calculate 

Ad> 


\Vemf\ 


At 


20.1  mV 


This  is  the  answer  to  part  (b).  In  some  sense,  this  is  the  average  emf  generated  in  the  coil 
over  the  200  ms  period  of  interest;  however,  more  precisely,  this  is  merely  the  best  estimate 
of  the  instantaneous  emf  generated  during  that  time,  given  the  limited  information  about 
the  time  dependence  of  the  magnetic  field  over  that  time. 
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[m0055]  [2] 

8.3-2 


Faraday’s  Law  says  that  the  emf  induced  in  a  this  loop  is 

iw  =  -!*(*> 

(since  a  loop  has  N  =  1  turn)  where 


<f>  =  /  B  •  ds 


Let  us  define  Vemf  to  be  across  the  resistor,  with  the  “+”  terminal  on  the  right  side  and  the 
”  terminal  on  the  left  side])]  Then: 


where  y o  is  the  location  of  the  left  side  of  the  loop.  ThenH 


r+L/2 

ryo+w 

$  =  B0 

/  dx 

/  eaydy 

J  x=—L/2 

o 

II 

1] 


Next  we’re  going  to  want  to  take  the  time  derivative  of  <f>.  ffowever,  to  do  that  properly 
we  need  to  make  sure  we  identify  everything  in  the  above  expression  for  $  that  has  a  time 
dependence.  Only  y0  depends  on  time.  To  make  this  clear,  let  us  write  yo  —  ut  +  b,  which 
places  the  left  side  of  the  loop  at  y  =  b  at  time  t  —  0.  Now  we  may  write: 


m 


TD  T 

aut  ab  r  aiu 
o  C-  1C 

a 


1] 


Vemf(t)  =  ~jtm  =  B0Lueauteab  [1  -  eaw ] 


The  problem  statement  asks  us  to  assess  the  situation  when  the  left  side  of  the  loop  is  at 
y  =  0.5  m,  so  we  choose  t  =  0  and  b  =  0.5  m  (since  we  said  earlier  y0  —  ut  +  b).  Also  from 
the  problem  statement,  u  =  —250  m/s.  Thus: 


Vemf(t  =  0)  =  =  0)  =  B0Lueab  [1  -  eaw]  =  -7.60  A 

Finally,  the  current  is  simply  this  divided  by  R,  which  is  —3.04  A.  Because  we  chose  the  “+” 
terminal  to  be  on  the  right,  the  reference  direction  for  current  must  be  counter-clockwise 
(i.e.,  this  is  the  necessary  direction  for  positive  current  to  dissipate  positive  power  in  the 

^^Not  the  only  way  to  do  it!  This  choice  is  arbitrary.  Choosing  the  opposite  reference  polarity  should  give 
you  the  exact  same  answer  as  long  as  you  follow  through  correctly. 

2 Once  you  complete  this  integration,  it’s  a  good  time  for  a  units  check! 

3Once  you  complete  this  differentiation,  it’s  a  good  time  for  a  units  check! 
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resistor,  or  alternatively  you  may  think  of  this  as  being  the  necessary  reference  direction  for 
the  loop  to  behave  as  a  power  source).  We  have  found  that  the  current  is  negative  with 


respect  to  this  reference  direction;  therefore,  the  induced  current  is  3.04  A,  clockwise 


You  can  check  to  make  sure  you  got  the  correct  current  direction  by  using  Lenz’s  Law.  If 
the  current  is  flowing  clockwise,  then  the  induced  magnetic  field  in  the  loop  is  in  the  —  z 
direction.  The  impressed  magnetic  flux  is  increasing,  since  the  loop  area  is  constant  and 
the  magnetic  field  in  the  loop  increases  in  the  +z  direction  as  the  loop  slides  to  the  left. 
Therefore,  the  induced  current  is  acting  to  oppose  the  change  in  the  impressed  magnetic 
flux,  as  is  required  by  Lenz’s  Law. 
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[m0055]  [3] 

8.3-3 


(a)  Let’s  stick  with  the  SI  system  of  units.  Since  “1”  in  the  expression  for  B  appears 
to  be  unitless,  Bq  must  have  units  of  T  (or  Wb/m2,  if  you  prefer)  .  Since  t  is  in  s  and  “1” 


appears  to  be  unitless,  k  must  have  units  of  1/s  .  Since  a  has  units  of  m  and  equals  vt,  v 
must  have  units  of 


m/s 


(b)  Here’s  Faraday’s  Law: 

v,{t)  =  -Nji  (8.1) 

where  N  —  1  since  it  is  a  loop  (not  a  coil),  and  the  magnetic  flux  is 

<F  =  /  B  •  ds  =  f  [z B0  (1  +  kt)]  ■  [— z ds]  =  —B0  (1  +  kt)  /  ds  (8.2) 

Js  Js  Js 

Note  that  ds  is  in  the  — z  direction.  This  is  from  the  right-hand  rule  (of  calculus),  in  which 
your  thump  is  along  the  loop  (not  the  gap)  pointing  along  the  direction  from  the  ”  terminal 
to  the  “+”  terminal.  Since  S  represents  the  surface  defined  by  the  loop  (actually  any  surface 
defined  by  the  loop,  but  we’ll  keep  it  simple...),  we  have 


<F  =  —  B0  (1  +  kt)  (tt a2) 

Substituting  a  =  vt  and  expanding  into  two  terms  we  get 

(8.3) 

$  =  — B0TTV2t 2  —  B0knv2t 3 

So  Faraday’s  Law  says: 

(8.4) 

Vg(t)  =  2B0nv2t  +  3B0kTTV2t2 
which  is  more  compactly  written  as: 

(8.5) 

Vg{t)  =  Bq-kv2  (2 1  +  3 kt2) 

(8.6) 

Good  time  for  a  units  check:  Can  you  confirm  that  the  result  is  dimensionally  correct? 


(c)  The  first  problem  is  to  determine  the  motional  and  transformer  emf,  so  we  should  be 
clear  on  what  we  mean  by  these  terms.  Motional  emf  is  the  contribution  to  the  total  emf 
which  is  associated  with  changes  in  the  size,  shape,  or  orientation  of  the  surface  through 
which  the  magnetic  field  lines  are  linked.  Transformer  emf  is  the  contribution  to  the  total 
emf  which  is  associated  with  changes  in  the  magnetic  held. 

With  that  in  mind,  let’s  consider  an  incorrect  solution:  You  can’t  set  v  —  0  and  call  the 
result  (in  this  case,  zero)  the  transformer  emf.  This  is  wrong  because  if  v  were  equal  to  zero 
for  some  radius  a  >  0,  then  the  calculated  emf  would  be  potentially  non-zero  because  the 
magnetic  held  is  still  time-varying. 

However,  it  is  true  that  setting  k  —  0  yields  the  motional  emf  (you  can  verify  this  for 
yourself  after  reading  through  this  solution),  although  this  is  pretty  hard  to  justify,  especially 
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since  we  just  saw  that  arbitrarily  setting  constants  to  zero  is  dangerous.  So,  if  you  go  that 
approach  you  must  be  able  to  explain  why  this  is  reasonable. 

With  all  this  in  mind,  here  are  three  reasonable  ways  to  get  a  solution: 

1.  You  could  calculate  the  motional  cmf  and  transformer  emf  using  Faraday’s  Law,  but  in 
two  separate  steps;  that  is,  work  out  the  emf  for  a  static  loop  (to  get  the  transformer 
emf)  and  then  for  a  static  held  (to  get  the  motional  emf). 

2.  You  could  calculate  motional  cmf  from  magnetostatics  (i.e.,  assume  a  static  magnetic 
held),  then  subtract  this  result  from  your  answer  for  part  (b)  to  get  the  transformer 
cmf. 

3.  The  approach  followed  below,  which  is  probably  best  since  it  deals  directly  with  the 
concepts  of  time-varying  loop  size  vs.  time-varying  magnetic  held. 

Here  we  go: 

m  =  -It‘!>  =  -ItIsB'ds  =  -7tIs '  [~ids]  (8J) 

where  B(t)  =  B0(  1  +  kt );  i.e.,  the  scalar  component  of  the  magnetic  held.  Continuing: 

vg(t)  =  +J  B(t)  (  ds]  =  +J  W)A(t)]  (8.8) 

s 

where  A(t)  =  7ra2(f);  i.e.,  the  area  of  the  loop.  The  value  in  setting  the  problem  up  this 
way  is  that  we  now  have  the  magnetic  held  and  the  loop  area  set  up  as  distinct,  identihable 
factors  in  the  solution.  Now  we  differentiate  using  the  chain  rule: 

va(t)  =  A(t)  +  B(t)  (8.9) 

Now  we  see  clearly  that  the  hrst  term  is  the  transformer  emf  and  the  second  term  is  the 
motional  emf.  Let’s  label  these  Vgr(t )  and  V™  respectively.  Now: 

K®  = 

Vgm(t)  =  Bit)  ^ A(t )  =  B(t)  ■  2na  ■  v  (8.11) 

y  I at 

Note  that  now  we  get  the  “expected”  result  when  we  set  k  =  0  and  then  v  —  0.  Also  note 
that  transformer  emf  depends  on  loop  area,  and  motional  emf  depends  on  loop  perimeter  - 
you  might  have  suspected  this  based  on  other  problems  you  have  encountered. 

OK,  now  we’re  ready  to  wrap  up.  To  hnd  out  when  the  contributions  of  the  transformer 
emf  and  motional  cmf  are  equal,  we  set  the  above  expressions  equal  and  solve  for  time  t  =  teq. 
Here  we  go: 

B0k  ■  it  a2  =  Biteq)  ■  2ira  ■  v  (8-12) 

Noting  that  a  =  vt  and  B(t)  =  B0  (1  +  kt)\ 

B0k  ■  iTV2tlq  =  B0  (1  +  kteq)  ■  2i Tv2teq  (8.13) 


— Bit)  Ait)  =  B0k  ■  tt a2 


(8.10) 
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Solving  the  above  expression  we  find  that  teq  =  —  2/k  (good  time  for  a  units  check,  by  the 
way).  Negative  times  don’t  really  make  sense  in  the  problem  (how  can  the  loop  have  zero 
radius  and  before  that  negative  radius?),  so  the  transformer  emf  and  motional  emf  are  never 
equal  for  k  >  0  (i.e.,  magnetic  field  magnitude  increasing).  However,  teq  =  +2/\k\  when 
k  <  0.  So,  the  answer  to  the  problem  is: 

at  t  =  —2/k,  and  only  if  the  magnitude  of  the  magnetic  field  is  decreasing. 
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[m0056]  [1] 

8.4-1 


The  measured  voltage  is  the  “transformer  emf”  induced  by  the  magnetic  flux  through  the 
loop.  According  to  Faraday’s  Law,  we  have  in  general  that: 

W  =  -*£*  = -"4  jf  b-* 

Here,  the  number  of  turns,  TV,  is  1;  ds  is  perpendicular  to  the  loop  in  the  direction  determined 
by  the  Stoke’s  Law  convention  (i.e.,  according  to  the  reference  polarity  chosen  for  Hprnf), 
and  B  can  be  written  as 

B  =  b  B(t)  =  bi?o  sin(27r/t  +  a ) 

where  b  is  simply  a  unit  vector  indicating  the  direction  of  B.  We  know  b  is  a  constant  with 
respect  to  position  because  B  was  specified  to  be  a  uniform  magnetic  field.  Thus,  we  have 
for  any  particular  orientation  of  the  loop  ds: 


V, 


emf 


d 

A  f 

B0 

—  sin(27r  ft  +  a) 

h  ■  ds 

|  Vgmfl  is  maximized  when  b  and  ds  point  either  in  the  same  direction,  or  in  exactly  opposite 
directions.  In  this  case,  the  magnitude  of  the  quantity  in  the  rightmost  square  brackets  is 
simply  the  area  of  the  loop,  A,  which  here  is  0.0314  m2  according  to  the  problem  statement. 
For  this  orientation,  we  have: 


I  ^emf  I  =  B o  [27t /  cos(27r/t  +  a)\A  . 

The  above  quantity  is  maximized  when  cos(27 rft  +  a)  =  1,  which  corresponds  to  the  peak 
magnitude  of  I/emp  which  is  one-half  of  the  peak-to-peak  magnitude.  Thus: 


B02irfA  =  0.5  x  20  mV™,  ,  thus: 


Bn  = 


0.5  x  0.02  V. 


pp 


2vr(100  x  103  Hz) (0.0314  m2) 


=  507  nT 
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8.4-2 


The  induced  potential  is 

d 

V  =  -IV— 4>  (8.14) 

where  4>  is  the  magnetic  flux  through  the  loop.  Since  the  magnetic  field  is  spatially  uniform 
and  the  loop  (in  each  case)  is  fixed,  we  may  express  the  flux  in  the  following  simple  form: 


$  =  AB0  cos  (cut  +  0) 


(8.15) 


where  A  is  the  area  of  the  loop,  Bo  is  a  constant  having  units  of  Wb/m2  (resulting  from 
the  dot  product  of  the  magnetic  flux  density  and  the  normal  to  the  loop),  and  c o  and  0 
are  the  angular  frequency  and  phase,  respectively,  of  the  sinusoidally-varying  flux.  It  is  not 
necessary  to  know  B0l  u,  or  0,  as  we  shall  see  in  a  moment.  Returning  to  the  Equation  18. 141 
we  see 

V  —  +NAB0lu  sin  (ut  +  0)  (8.16) 

Therefore  the  peak  potential  is  simply  N ABqLO.  So  here’s  the  situation: 

y't2)  nvavbw 

l/(«  NWAWBau  1  ' 

pk 


where  the  superscripts  indicate  before  (i.e. ,  one  turn  circular  loop)  vs.  after  (i.e. ,  two-turn 


square  loop).  Thus,  =  15  V  and  we  seek  V00  Since  neither  B0  nor  uj  change  between 


(2) 


pk 


the  two  scenarios,  we  have: 


Solving  for  V^\  we  find: 


K 


(2) 


pk 


pk 


n^a^ 


15  V  NWAW 


(8.18) 


2  •  (0.2  m)2 
•  7T  (0.1  m)2 


38.2  V 


(8.19) 
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8.5-1 


In  the  original  scenario: 


V2  N2 
v !  “ PN 1 


where  p  —  ±1  depending  on  the  relative  orientation  of  the  windings  (p  —  —  1  for  the  example 
shown  in  the  book).  Let  Vf  be  the  new  potential  on  the  secondary  coil.  From  Faraday’s 
Law: 

V'  =  iv4*'2 

where  $2  is  the  magnetic  flux  through  the  secondary  coil  after  the  modification.  Note 

$2  —  — 

2  2 

since  the  secondary  coil  now  intersects  only  half  the  flux  it  did  previously.  Subsequently, 

=  PY 


Now: 


V2  =  p\N2^  =p\j^  i 


dt  7  r2  Nx 

From  the  problem  statement,  N\  =  200  and  N2  =  300.  Also  p  =  —  1  for  the  example  shown 
in  the  book.  Therefore, 


V2'=  --Li 


with  the  sign  depending  on  the  relative  orientations  of  the  coil  windings. 
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[m0031]  [2] 

8.5-2 


The  transformer  is  an  application  of  Faraday’s  law,  which  is  intrinsic  to  the  Maxwell-Faraday 
Equation: 

V  x  E  =  — 4b 
ot 
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8.7-1 


From  the  problem  statement: 

B0  =  2  T, 

A  =  7T  (2  m2)  =  47t  m2, 

the  peak  value  of  Vr  is  5  V,  and 

Vrit  —  0)  =  0  and  increasing. 

Since  the  loop  is  rotating  in  a  static  uniform  magnetic  field,  Vt  must  be  sinusoidally- varying. 
A  general  form  for  this  variation  is 


Vr(t)  =  AB0uj  cos  (ut  +  VO 


where  u  =  27t/  is  the  angular  frequency  of  rotation  and  V*  is  an  as-yet  unknown  phase  offset. 
However,  it  is  known  that 

Vr(t  —  0)  =  AB0u  cos(V0  =  0 

so  must  be  either  ir/2  or  3n/2.  Since  Vr  is  increasing  at  t  —  0,  V’  must  be  3n/2. 
Furthermore,  we  know  the  peak  value  of  Vt-  In  the  context  of  the  general  form,  we  find: 


Vr  =  ABqoo  =  5  V  at  maximum 


and  therefore,  c o  =  0.199  rad/s. 


Putting  this  all  together: 


VT{t)  =  (5  V)  cos  (  [0.199  rad/s]  t  + 
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8.7-2 


Assigning  symbols  to  quantities  identified  in  the  problem  statement:  B0  =  2  T,  A  =  4  m2, 
and  Vr\pk  =  5  V.  Recall: 

VT  =  2nf0AB0b  •  p(t)  (8.20) 

where  /o  is  the  frequency  of  rotation,  b  is  the  direction  of  the  magnetic  field,  and  p(t)  is  a 
unit  vector  that  lies  in  the  x  —  y  plane  and  in  the  plane  of  the  loop,  rotating  with  the  loop. 
Note  the  maximum  magnitude  of  b  •  p(t)  is  simply  1.  Thus: 


Vr,pk  —  27r/oA-E>o  —  5  V 


Solving  for  /0: 


/o 


5  V 

2tt  •  (4  m2)  •  (2  T) 


99.5  mHz 


(8.21) 

(8.22) 
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8.9-1 


From  the  problem  statement,  the  electric  held  intensity  is 

E  =  y  (3  V  rif1  s“2)  t2 

(The  fact  that  this  is  electric  field  intensity  can  be  confirmed  using  dimensional  analysis.) 
The  displacement  current  density  in  free  space  is 


^  =  ye0  (6  V  m  1  s  2)  t  =  y  (53.1  pC  m  2  s  2)  t  = 


dt 


y  53.1 


pA 

m2  s 
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Chapter  9 


Plane  Wave  Propagation  in  Lossless 
Media 
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9.1-1 


The  general,  time-domain,  differential  form  of  Ampere’s  Law  is: 


<9D 

V  x  H  =  J  +  — — 
at 

The  relationship  between  these  quantities  and  the  phasor  representation  of  the  same  quan¬ 
tities  is: 

H  =  Re  \  Uejut 


J  =  Re|je^|  ,  and 
D  =  Re  { i)ejut 


Now  substituting  these  quantities  into  Ampere’s  Law  we  have: 


V  x 


Re  \  Hejujt 


=  Re  \  J  ejujt  \  + 


Re  \  BejuJt 


The  order  of  the  “Re”  operator  and  any  linear  real-valued  operator  can  be  exchanged  (see 
the  textbook  section  on  phasors  for  a  proof  of  this).  Taking  advantage  of  this  in  the  first 
and  last  terms,  we  obtain: 


Re(Vx 


Hej0Jt 


=  Re  l  Jejut 


+  Re  <!  ^  [Ddwt 


Note  that  the  curl  (“Vx”)  operator  operates  only  on  position,  and  not  on  time.  Thus,  we 
may  rewrite  the  first  term  as  shown  below: 


Re 


V  x  H 


jwt  \  =  Re  \  3ejujt 


+Re<m 


Dejut 


Note  also  that  partial  derivative  in  the  last  term  operates  only  on  time,  whereas  D,  being  a 
phasor,  is  independent  of  time.  Therefore,  the  partial  derivative  operates  only  on  the  factor 
eJujt ,  and  we  have: 


Comparing  terms 
is: 


Re 


V  x  H 


s^} 


above,  we  find  that 


=  Re  jje^j  +  Re  {D  •  juejut}  . 

the  phasor  expression  of  Ampere’s  Law  that  we  seek 


V  x  H  =  J  +  ju D  . 
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9.1-2 


The  expression  worked  out  in  the  book  is 

V  x  E  =  -juB 

B  is  a  flux  density,  so  we  use  B  =  pH  to  obtain: 

V  x  E  =  —  jcu/iH 
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9.2-1 


Here  are  Maxwell’s  Equations  for  source-free  regions  in  terms  of  E  and  H  (only)  in  dif¬ 
ferential  form: 

V  ■  E  =  0  (9.1) 

V  x  E  =  — jcu/iH  (9.2) 

V  -  H  =  0  (9.3) 

V  x  H  =  +juje E  (9.4) 

The  equation  we  seek  must  yield  solutions  for  H  which  satisfy  at  least  the  last  three  of  the 
above  four  equations.  We  begin  by  taking  the  curl  of  Equation  19.41 

V  x  x  Hj  =  V  x  (^-jcoeE^j  =  +jcue  (V  x  E^  (9.5) 

On  the  right,  we  can  substitute  for  VxE  using  Equation  19.21 

+  jcoe  x  E^  =  +joje  =  +o ;2/reH  (9.6) 

On  the  left,  we  invoke  the  vector  identity 

VxVxA  =  V(V-A)—  V2A  (9.7) 

to  obtain 

VxVxH  =  v(V'H)-  V2H  =  -V'2H  (9.8) 

where  we  have  used  Equation  19.31  to  eliminate  the  V  •  H  term.  Substituting  back  into 
Equation  19.51  and  rearranging  terms  we  have 

V2H  +  cu2/ieH  =  0  (9.9) 

Now  substituting  f3  =  ojy/JIe\ 

V2  H  +  /32H  =  0  (9.10) 

This  is  the  homogeneous  wave  equation  for  H. 
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9.2-2 


The  wave  equation  for  E  is 


V2E  +  /32E  =  0 


Also,  we  know  that 


_  9°°  _  tt/2  rad  _ 
1  m  1  m 

which  may  also  be  expressed  simply  as  7t/2  m  .  So 


7 r 
2 


rad/m 


V2E  +  (2.467  m”2)  E  =  0 
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9.4-1 


(a)  The  wave  equation  for  E  is  V2E  +  /32E  =  0.  In  cylindrical  coordinates,  E  =  pEp  + 
(fiE'p  +  z Ez.  The  Laplacian  operator  in  cylindrical  coordinates  is: 

2  _  {  d_\  1  d2  d2 

pdp  \  dp )  p2  dcj)2  dz 2 

Thus,  we  have  for  the  three  components  of  E: 


1  d 

pdp 

/  q  _ \  ]_  3^  _  3^  ~  ~ 

\pdpE'‘)  +  p--arEp  +  dziEp  +  p  E°  - 0 

1  d  / 

pdp  \ 

fdpE*)  +  p--  a +  dz2E*  +  P  E*  ~ 0 

^  |c§- 

T - 1  1 

(  d  ~  \  l  d2  ~  d2  ~  ~ 

\dpEz)  +  p2  d^Ez  +  dz2Ez  +  ^  Ez  ~  ° 

(b)  If  E  has  no  component  in  the  p  or  0  direction,  then  Ep  =  E $ 
one  equation: 


1  d_ 

pdp 


1  d2  ~  d2  ~  ~ 

H - 2 

pz  dcp  dzz 


0  and  we  are  down  to 

0 


If  E  is  uniform  in  0  and  z,  then  dEz/d(f)  =  dEz/dz  =  0,  so  the  second  and  third  terms  in 
the  above  equation  are  zero.  This  leaves  us  with: 


1  d_ 

pdp 


+  P2ez 


0 
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9.4-2 


First  note  that  xcos  0  +  y  sin  0  is  a  unit  vector;  i.e. , 


|xcos  0  +  y  sin  0|  —  cos2  0  +  sin2  0  =  1  (9.11) 

Next,  recall  that  E  x  H  is  in  the  direction  of  propagation,  which  in  this  problem  is  +z. 
Therefore  the  direction  of  H  is  z  x  E.  (If  this  is  not  clear,  think  of  E,  H,  and  direction  of 
propagation  forming  a  cartesian  coordinate  system  with  E  analogous  to  x,  H  analogous  to 
y,  and  direction  of  propagation  analogous  to  z.)  Thus  the  direction  of  H  is 

z  x  (x  cos  0  +  y  sin  0)  =  —  x  sin  0  +  y  cos  0  (9-12) 

and  this  is  a  unit  vector  since  xcos0  +  ysin0  is  a  unit  vector  which  is  perpendicular  to  z. 


The  magnitude  of  H  is  |E|  / 77,  where  r]  in  this  case  is  a/puT^o  —  377  fl.  Therefore  the  mag¬ 
nitude  of  H  is  (2  /iV / m)  /rj0  =  5.31  nA/m. 

Putting  this  all  together,  the  magnitude  and  direction  of  the  associated  magnetic  held  is: 


(—x sin  0  +  y  cos  0)  5.31  nA/m 


(9.13) 
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9.5-1 


(a)  There  are  several  ways  to  figure  this  out.  One  way  is  to  start  in  phasor  represenation, 
in  which  the  field  is: 

eE0ej{ax+by+cz) 

where  a  =  1  rad/m,  b  =  2  rad/m,  and  c  =  3  rad/m.  Note: 

eE0ej{ax+by+cz)  =  eE0e+jaxe+jbye+jcz 


To  see  what’s  going  on  here,  consider  a  simpler  version  of  the  above  equation,  where  the  last 
two  terms  are  omitted: 

eE0e+jax 

In  this  case,  the  answer  we  are  looking  for  would  be  k  =  —  x,  and  in  fact  we  also  see  that 
the  wavenumber  [3  —  a;  i.e. ,  1  rad/m.  Similarly,  if  we  had  just 

eE0e+jby 


then  the  answer  would  be  k  =  —  y,  with  /3  —  b  —  2  rad/m.  From  this  we  can  infer  that  a 
vector  (not  necessarily  a  unit  vector)  that  points  in  the  direction  of  propagation  in  this  case 
is: 

k  =  —ax  —  by  —  cz 

Thus,  the  corresponding  unit  vector  is 


Giving: 


k  —ax  —  by  —  cz 

kf  =  V(~a)2  +  ( -b )2  +  (-c)2 


k 


— 0.267x  —  0.534y  —  0.802z 


(b)  In  the  above  analysis,  we  find  |k|  =  3.74  rad/m.  Following  the  reasoning  above,  this  is 
simply  the  wavenumber  (3.  Thus,  the  wavelength  is: 


A 


2v r 

J 


1.68  m 


(c)  Since  this  is  free  space,  and  since  we  know  the  phase  velocity  in  free  space  is  c  = 
3.0  x  108  m/s,  we  also  know  the  frequency,  which  is  c/A= 


179  MHz 
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9.5-2 


From  the  problem  statement,  the  direction  of  propagation  k  =  —  x  and  H  points  in  the 
+y  direction.  From  the  plane  wave  relationships: 


Therefore,  E  points  in  the  - 

E  =  —77k  x  H 

(— x)  x  y  =  +z  direction. 

141 


[m0039]  [3] 

9.5-3 


From  the  problem  statement,  |E|  =  3  V/m  and  er  =  2.  Since  plastics  are  non-magnetic: 


V  = 


V’O  Vo 


ere0 


where  r]0  =  376.7  fl.  Therefore, 


I  HI 


|E|  _  |E| 
V  Vo 


6r  — 


11.3  niA/m 
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9.7-1 


The  spatial  power  density  is 


S avp. 


3  W 
1  mm2 


3  W 
10-6  m2 


=  3  MW/m2 


Since  Save  =  |E|2  /2p,  and  since  rj  —  rjo  —  376.7  17  in  free  space: 


|E|  =  sj2r]0Save  =  47.5  kV/m 
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